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results) 
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normal-force coefficient, -- 

co^l y 

local normal-force coefficient per unit length 


pressure coefficient. 


side-force coefficient, v 

Q oo ^ y 

body cross-section diameter 


cross force per unit length along body length 
potential cross force per unit length along body length 
viscous cross force per unit length along body length 
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longitudinal and transverse apparent mass coefficients 
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free-stream Reynolds number, 
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body volume 
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P 
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reference length (generally taken as d for the comparisons of computed with 
experimental results) 

axial distance from body nose 

axial distance from body nose to aerodynamic normal-force center (center of 
pressure) 

axial distance from body nose to centroid of body planform area 

axial distance from body nose to pitching-moment reference center 

axial distance from body nose to aerodynamic side-force center 

angle of attack 

angle of sideslip 

wing planform semiapex angle 

crossflow drag proportionality factor (ratio of crossflow drag for a finite-length 
cylinder to that for an infinite-length cylinder) 

viscosity coefficient of air 

density of air 

angle of roll about body longitudinal axis 
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PREDICTION OF STATIC AERODYNAMIC CHARACTERISTICS 
FOR SLENDER BODIES ALONE AND WITH LIFTING 
SURFACES TO VERY HIGH ANGLES OF ATTACK* 

Leland Howard Jorgensen 

Ames Research Center 


SUMMARY 

An engineering-type method is presented for computing normal-force and pitching- 
moment coefficients for slender bodies of circular and noncircular cross section alone and 
with lifting surfaces. In this method, a semiempirical term representing viscous-separation 
crossflow is added to a term representing potential-theory crossflow. In computing Cry and 
r for bodies alone, slender-body theory is used for the term representing the potential 
crossflow. For bodies with thin wings and tails, the linearized potential method of Nielsen. 
Kaattari, and Pitts, modified for high angles of attack, is used. 

For many bodies of revolution, computed aerodynamic characteristics are shown to 
agree with measured results for investigated free-stream Mach numbers from 0.6 to 2.9. The 
angles of attack extend from 0° to 180° for M ^ = 2.9 and from 0° to 60° for = 0.6 
to 2.0. 

For several bodies of elliptic cross section, measured results are also predicted reason- 
ably well over the investigated Mach number range from 0.6 to 2.0 and at angles of attack 
from 0° to 60°. As for the bodies of revolution, the predictions are best for supersonic Mach 
numbers. 

For body-wing and body-wing-tail configurations with wings of aspect ratios 3 and 4, 
measured normal-force coefficients and centers are predicted reasonably well at the upper 
test Mach number of 2.0. However, with a decrease in Mach number to M a 0 - 0.6, the 
agreement for Cat rapidly deteriorates, although the normal-force centers remain in close 
agreement. 

For M 00 = 0.6, 0.9, and 2.0 and angles of attack of 10°, 20°, 30°. 40°, and 50°. 
vapor-screen and oil-flow pictures are shown for many body, body-wing, and body-wing-tail 
configurations. When separation and vortex patterns are asymmetric, undesirable side forces 
are measured for the models even at zero sideslip angle. 

These side forces can be significantly affected by changes in Mach number, nose 
fineness ratio, nose bluntness, and body cross section. Generally, the side-force coefficients 
decrease or vanish with the following: increase in Mach number, decrease in nose fineness 
ratio, change from sharp to blunt nose, and flattening of body cross section (particularly the 
body nose). Additions of afterbody strakes. wings, or wings plus tail produce much smaller 
or no appreciable effects. 


♦Formerly issued as NASA TM X-73,123 with limited distribution. Reissued for unlimited distribution. 



CHAPTER 1 
INTRODUCTION 


Over the last several years, high angle-of-attack aerodynamics has increased in impor- 
tance because of the demand for greater maneuverability of space shuttle vehicles, missiles, 
and military aircraft (both manned and remotely piloted). Until recently there has been a 
general lack of analytical methods and aerodynamic data suitable for use in the preliminary 
design of most advanced configurations for flight to high angles of attack over a wide range 
of Mach and Reynolds numbers. There has been, however, considerable research leading to 
the development of methods for predicting the static aerodynamic characteristics of simple 
shapes, primarily slender bodies of revolution. 

Prior to the work of Allen in 1949 (ref. 1), most analytical procedures for computing 
the aerodynamic characteristics of bodies were based on potential-flow theory, and their 
usefulness was limited to very low angles of attack. Allen proposed a method for predicting 
the static longitudinal forces and moments for bodies of revolution inclined to angles of 
attack considerably higher than those for which theories based only on potential-flow con- 
cepts are known to apply. In this method, a crossflow lift attributed to viscous-flow separa- 
tion is added to the lift predicted by potential theory. This method has been used quite suc- 
cessfully to compute the aerodynamic coefficients of inclined bodies (e.g., refs. 1-6), 
although most data available for study until 1961 were for bodies at angles of attack below 
about 20°, and the formulas were initially written to apply only over about this angle-of- 
attack range. 

In 1961, Allen’s concept was adapted by Jorgensen and Treon (ref. 7) for computing 
the normal-force, axial-force, and pitching-moment coefficients for a rocket booster 
throughout the angle-of-attack range from 0° to 180°. Reasonable agreement of theory with 
experiment was obtained for a test model of the rocket booster over the Mach number range 
from 0.6 to 4. The Allen concept was further applied by Saffell, Howard, and Brooks 
(ref. 8) in 1971 in a computer-programmed method for predicting the static longitudinal 
aerodynamic characteristics of low aspect-ratio missiles operating at angles up to 1 80°. 

In 1958, a method tor computing the aerodynamic characteristics for bodies of non- 
circular cross section was proposed by Jorgensen (ref. 6). In this method, normal-force and 

pitching-moment coefficients (C v and C m ) are computed by Allen’s formulas for the 

o n 
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equivalent body of revolution which has the same axial distribution of cross-sectional area as 
the noncircular body. Then the values of Cjy and C fn for the noncircular body are computed 
from Cy/ C.y and C m /C m ratios determined from apparent mass coefficients (i.e., from 
slender-body theory). Good agreement of theory with experiment (ref. 6) was obtained by 
this procedure for bodies of elliptic cross section at the conditions investigated (a/b = 1 
to 2,0 — 0° and 90 °,M 00 = 2 to 4, and a = 0° to 20°). 

Early in 1973, the Allen concept was again applied by Jorgensen (ref. 9) to further 
develop an engineering-type procedure for computing normal-force, axial-force, and 
pitching-moment coefficients for slender bodies of circular and noncircular cross sections at 
a = 0° to 180°. The Cy and C m formulas were written, however, for a body whose cross- 
sectional shape remains constant over the body length, but the cross-sectional area, ot 
course, is allowed to vary. 

In 1973, Jorgensen (refs. 10 and 1 1) rewrote the Cy and C fn expressions to apply for 
the general case of a slender body alone or with lifting surfaces where the cross-sectional 
shape, as well as the cross-sectional area, is allowed to vary along the body length. For the 
special case of winged-elliptic cones, simplified expressions for Cy and C... were also pre- 
sented. Good agreement between predicted and experimental results was shown. However, 
experimental results available for comparison with the method were limited to simple ellip- 
tic bodies and winged-elliptic cones at angles of attack less than about 20° and Mach 
numbers only from 2 to 4. Thus, it was concluded that additional comparisons for these and 
more diverse configurations were needed at higher angles of attack and lower Mach numbers 
to determine validity limits for the method. 

In 1972-73 there was a great need to enlarge the relatively small data base for bodies 
alone and in combination with wings and tails at subsonic, transonic, and supersonic Mach 
numbers. This need still exists today, but it has been alleviated significantly by some recent 
introductory investigations into this high a field (e.g., refs. 12—21). Most studies have been 
initiated primarily for application to missile aerodynamics. Two recent studies, however, 
have been more aircraft oriented in that they have been directed toward the determination 
of experimental aerodynamic characteristics for slender bodies with thin wings (ret. 20) and 
with wings and a tail (ref. 21). There is still need for study and analysis of much of this 
recently obtained data both for slender bodies alone and for bodies with wings at very high 


3 



angles of attack. As previously mentioned, the Jorgensen method should be further com- 
pared with experiment. Also, there is an urgent need for continued study of unwanted side 
forces and yawing moments which have been shown to develop for models at high a values 
with (3 = 0°. These unwanted side forces and yawing moments, which probably lead to 
undesirable stall/spin characteristics, have been measured on noses alone, bodies alone, 
bodies with wings, and bodies with wings and a tail (e.g., refs. 12-21). 

In view of the foregoing, the present study was initiated to accomplish the following 
objectives: 

1. Review and extend the derivation of Jorgensen’s engineering-type method 
(refs. 9-11) for computing the normal-force and pitching-moment characteristics of slender 
circular and noncircular bodies alone and with thin lifting surfaces. 

2. Assess the method for predicting the aerodynamic characteristics of missile-like 
bodies of revolution. 

3. Assess the method for predicting the aerodynamic characteristics of bodies of ellip- 
tic cross section. 

4. Assess the method for predicting the aerodynamic characteristics of bodies with 
thin wings alone and in combination with a tail. 

5. Present and discuss visual observations of the vortex flow fields over models of 
bodies alone and in combination with wings. 

6. Present and discuss data pertaining to the origin and alleviation of undesirable side 
forces and yawing moments associated with high a flight at (3 = 0°. 

In connection with objectives (5) and (6), the present study is limited in scope to an 
experimental study of the vortex flow fields and side forces associated with bodies and 
body-wing combinations. It is firmly believed that a larger experimental foundation must be 
acquired at this time to aid in the further development of flow-field modeling techniques 
that can be applied with confidence in preliminary design studies. This is not to say that 
important semiempirical modeling techniques have not already been studied. Most of these 
techniques have been based on an impulsive flow analogy and applied to slender bodies of 
revolution. In this analogy, it is assumed that the crossflow plane is swept uniformly down 
the length of an inclined body at the rate cos a. Then there is an analogy made between 
the leeward flow field and the developing wake behind an impulsively started cylinder in 
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crossflow. This analogy was initially suggested by Allen and Perkins (ref. 2) and applied by 
Kelly (ref. 22). 

The impulsive flow analogy probably has been used most elegantly in recent studies by 
Wardlaw (refs. 23 and 24). In one of his latest studies (ref. 24), he simulated the viscous 
crossflow plane by superimposing a large number of point vortices on the potential solution 
for flow about a cylinder. For some pointed-nose bodies, he has shown qualitative agree- 
ment with experimental observations, but even the latest Wardlaw method still must be 
considered to be in an early exploratory phase. Others have used the impulsive-flow analogy 
and attempted to model the crossflow field (e.g., refs. 25-32). 

Attempts also have been made to solve the crossflow field problem with the Navier- 
Stokes equations (refs. 33 and 34), but these studies have been applied only to bodies at 
supersonic speeds and moderate to low angles of attack. Computer times appear to become 
prohibitive for most practical cases, and turbulence modeling is not yet demonstrated. 

Because of many deficiencies in the understanding of the physical flow fields around 
bodies and wings at very high angles of attack, further experimental investigation is required 
before much reliability can be realized from further analytic modeling. 
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CHAPTER 2 

DERIVATION OF BASIC METHOD FOR COMPUTING Cy AND C m 

CHARACTERISTICS 

Here we first review the derivation of some basic equations for computing the 
normal-force and pitching-moment coefficients, Cj y and C m , for slender bodies of revolu- 
tion to very high angles of attack (section 2.1). Then we extend the derivation to obtain 
more general Cy and C m equations that can be used for slender bodies of circular and 
noncircular cross section alone and with thin lifting surfaces (section 2.2). Necessary empir- 
ical input values ot crossflow drag coefficient for circular cylinders are presented in sec- 
tion 2.3, and theoretical formulas to obtain input local normal-force coefficients for non- 
circular configurations are presented in section 2.4. Finally, in section 2.5, we briefly discuss 
the relative influence of derived potential-flow and viscous crossflow terms used in the basic 
equations to compute both Cy and C m as a function of angle of attack. 

2. 1 Bodies of Revolution 

In 1949, H. J. Allen (ref. 1) proposed a heuristic concept for predicting the static longi- 
tudinal forces and moments for bodies of revolution inclined to angles of attack consider- 
ably higher than those tor which theories based only on potential-flow concepts are known 
to apply. In this concept, a crossflow lift attributed to viscous crossflow separation is added 
to the crossflow lift predicted by potential-flow theory. For the potential-flow lift. Allen 
used the slender-body equation derived in 1923 by Max Munk for airship hulls (ref. 35). 

From the momentum consideration, Munk (ref. 35) showed that, for slender (high- 
fineness-ratio) bodies, the potential-flow cross force per unit length f at any station along 
the body is given by 

fp= (K 2 - K\ )q 00 sin 2a (2.1 ) 

where A is the cross-sectional area of the body at any axial distance x from the nose apex; a 
is the angle of attack; and A 2 and are, respectively, the transverse and longitudinal 
apparent mass coefficients tor the body. From calculations of A 2 and K ( made initially by 
H. Famb (ref. 36) for ellipsoids of various fineness ratios, Munk (ref. 35) has shown that 
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K 2 - K t is approximately unity for high fineness ratios. It has been customary to assume a 
value of unity for the factor K 2 - Ki, and it is assumed in this formulation. 

G. N. Ward (ref. 37) has shown that the section cross force f at small angles of attack 
acts at an angle midway between the normal to the axis of revolution of the body and the 
normal to the free-stream velocity (i.e., at an angle a/2). With this consideration, equa- 
tion (2.1 ) is multiplied by cos (a/2) to determine the normal-force distribution. At high a, the 
potential lift becomes small in comparison with the viscous crossflow lift, and it is of little 
practical consequence whether this low a multiplier [cos (a/2)] is retained or replaced with 
unity. 

Now consider the derivation of the crossflow lift attributed to viscous crossflow separa- 
tion. Allen (ref. 1), for a body of revolution of high fineness ratio, first treated each circular 
cross section as an element of an infinitely long circular cylinder of the same cross-sectional 
area. With this assumption, the local cross force per unit length due to viscosity is given by 

P V n 

f, - 2 ' c d n -f- < 2 - 2 > 

where r is the body radius at any station x from the nose apex; V n is the velocity normal to 
the longitudinal axis; p is the mass density; and is the steady-state crossflow (or normal) 
drag coefficient based on q„, diameter, and unit length at station*. Cj is a function of 
both the Mach number and Reynolds number components normal to the cylinder longitu- 
dinal axis. Hence, for a body at angle of attack, Cj is a function of 

u n 

V n 

M n =—= sin a (2.3) 


and 


Re n = Re sin a 


where 

v n = V oo sin « 


(2.4) 


(2.5) 


Here M^, Re, and are free-stream values, respectively, of Mach number, Reynolds 

number, and velocity; the speed of sound is denoted by a. 
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For this study, C i is taken as the steady-state crossflow drag coefficient for a circular 
u n 

cylinder. Other researchers (e.g., refs. 22-32) have assumed that the development of the 
crossflow with distance along an inclined body of uniform diameter is analogous to the 
development with time of the flow on a cylinder impulsively set in motion from rest. In this 
analogy, it is assumed that the crossflow plane is swept uniformly down the length of an 

inclined body at the rate cos a. When this analogy is assumed, the value of Cj starts 

a n 

from zero (at zero time), then increases to a maximum value about 25 percent higher than 
the steady-state value, then decreases to the steady-state value for laminar flow (see, e.g., 
sketch (a) from ref. 27). For turbulent flow, experiments conducted by Sarpkaya (ref. 27) 
indicate that the steady-state condition is reached almost at the start of the motion. Thus, 
for turbulent flow, both the steady-state assumption and the impulsive-flow assumption give 
the same result. However, for laminar flow there is a difference, but this difference is not 
studied in this report. 



Sketch (a) 

Because of spillage flow around the ends of a finite length cylinder, the value of Q is 

a n 

less than that for an infinitely long (truly two-dimensional) cylinder in the same free 
stream. Equation (2.2) should be multiplied by a proportionality factor rj, which is the ratio 
of C d n for a finite len gth cylinder to that for an infinite length cylinder. This factor rj, 
which approaches unity as the cylinder length to diameter approaches infinity, is given from 
experimental results (discussed later). In the practical use of tj, it is assumed that the £ jd of 
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the cylinder is the same as that for the body of revolution being considered. It is also 
assumed that the net-force effect of the front and rear end flows is approximately the same 
for both configurations. Thus, equation (2.2) is modified to give 


where 


f v = 2 wC d ^ sin 2 a 


= I p v 2 = - pv 2 — ? — 

2 P 00 2 P « sin 2 a 


sin 2 a 


( 2 . 6 ) 


(2.7) 


Following Allen (ref. 1), the potential-flow solution of Munk and the viscous crossflow 
solution are combined to determine the cross-force distribution along the body of revolu- 
tion. With the potential cross force per unit length acting at a/2 from the normal to the free- 
stream direction and the viscous cross force per unit length acting normal to the longitudinal 
axis, the total cross force per unit length normal to the longitudinal axis is given by 

f = f p cos~ + f v (2.8) 


= sin 2a cos ^ + 2r?Cy (sin 2 a)r (2.9) 

7oo 2d x a n 


From equation (2.9), the equations for the normal-force and pitching-moment coeffi- 
cients for bodies of revolution can be easily derived. With the normal force given by 

F n = f dx and the normal-force coefficient defined by Cjq - (F^q^y), we obtain 

o 


C N ~ 


sin 2a cos 




dA 

dx 


2ij Ctf sin 2 a /•£ 


r dx 


( 2 . 10 ) 


Likewise, with the pitching-moment coefficient defined by 


„ _ pitching moment _ 

m qZaJc 


f 

•V) 


f(x m - x)dx 




we obtain 


C, 


sin 2a cos (a/2) 


m 


A r X 


r 


CM 2vC <ln 

d? <x m ~ * )d * + — ipr 


sin 2 a 


/' 


r(x 


m 


x)dx 

( 2 . 11 ) 


9 



where A r is the reference area; X is the reference length; C is the body length; and x m is the 
axial distance from the nose to the pitching-moment reference center. 

General integrated expressions for bodies of revolution at a = 0° to 1 80° can be 
obtained from equations (2. 10) and (2.1 1). For the sign convention in sketch (b), Jorgensen 
(ref. 9) has written the following equations for the normal-force and pitching-moment coef- 
ficients; 

Cjy - sin 2a' cos sin 2 a' ; 0°<a<180° (2.12) 

Ay L* Ay 


r = 
^ m 


i 


V- A b (Z - X m ) 


A r X 


Si „ 2a' cos “ + „C dn 


sin* a 


0° < a < 90° (2.13) 

and 



90° < a < 180° (2.14) 

where A b is the body base area; A p is the planform area; V is the body volume; x ( . is the 
distance from the nose apex to the centroid of planform area; and a’ = a for 0°< a < 90° 
and a = 180°- a for 90° <a< 180°. 
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The axial distance from the nose apex to the aerodynamic force center is then given by 



When lift and drag coefficients are desired, they may be obtained at all a values from 
the conversion expressions: 

Cl = C/v cos a " C4 sin a (2.16) 

Cjy = Cjy sin a + C4 cos a (2. 1 7) 

Generally, it can be shown that for values of a well removed from a = 0° and 180° precise 

prediction of is not necessary to obtain reasonably accurate values of and Cjy 

Jorgensen (ref. 9) has suggested that, for rough engineering estimates, the axial-force 

coefficients for slender bodies can be approximated by 

Ca*Ca 0 cos 2 a ; 0° < a < 90° (2.18) 

^ ^ 0 

and 

Ca^Ca o cos 2 a ; 90°<a<180° (2.19) 

A ^a=i80 

Here cos 2 a is merely the ratio of the dynamic pressure in the axial direction to the dynamic 
pressure in the free-stream direction. A more precise method is also outlined in reference 9, 
along with procedures for determining C 4 _ q0 and C A _ I80 <> for often used conical and 
tangent ogive noses and flat bases. 

2.2 Bodies of Circular and Noncircular Cross Section 
Alone and With Lifting Surfaces 

For the general case of a body alone or with lifting surfaces where the cross-sectional 
shape can vary along the body length, procedures similar to those used for a body of 
revolution are assumed. Both the potential and viscous separation crossflow terms in equa- 
tions (2. 1 0) and (2. 1 1 ) for Cjy and C m are generalized further. 

The potential crossflow term is generalized by multiplying the value inside the integral 

by the ratio (C„/C„ ) - the ratio of the local normal-force coefficient per unit length C n 

0 SB 

for the desired cross-sectional shape to the similar coefficient C„ lor the equivalent circular 

n o 



shape having the same cross-sectional area. The necessary ratios can be determined from 
apparent mass coefficients (slender-body theory) for many cross-sectional shapes. In 1958, 
this procedure was first shown (ref. 6) to have merit in predicting experimental Cyy and C m 
results for slender bodies of various cross sections at supersonic speeds and angles of attack 
up to about 20°. However, because of very limited experimental data, further evaluation 
was not possible at that time. 


The viscous crossflow terms in equations (2.10) and (2.1 1) are further generalized by 

multiplying the values within the integrals by (C„/C„ ) . Here the local ratio of C„ to 

„ ° Newt n 

C n at each x station is assumed to be given by Newtonian impact theory. Cj remains as 
o a « 

the crossflow drag coefficient for the equivalent circular cylinder section. 

With equations (2.10) and (2.11) generalized as discussed, and for positive cL4/dx 


n _ sin 2a cos 

c N 


sin 2a cos (a/2) 


2 r?C^ sin 2 a 


x cos (a/2) C t^n 

N l IS 

- sin 2 a \ 

W (4 

' ° 'SR 


dA , ,, 

to {X m - x)dx 


r(x m - x)dx 


( 2 . 20 ) 


( 2 . 21 ) 


Newt 


In equations (2.20) and (2.21), the first terms (from slender-body theory) are not applicable 
as written for winged-body sections where the cL4/dx values are zero or negative, and 
procedures similar to those suggested in reference 38 probably should be used. Further 
adaptation of this method for use with body-wing and body-wing-tail configurations is 
considered in chapter 5. 

In the second term of equations (2.20) and (2.21), there is some experimental justifica- 
tion for formulating the ratio (C„/C„ ) from Newtonian theory and multiplying it by the 
available experimental or theoretical crossflow drag coefficient for the equivalent circular 
cross section C d For subcritical crossflow Mach and Reynolds numbers, Jorgensen 
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(ref. 11) has shown that CJC„ values from Newtonian theory agree reasonably well (but 

n ri Q 

somewhat fortuitously) with those from two-dimensional tests (refs. 39-43) of elliptic 
cross sections and square cross sections with rounded corners. Jorgensen’s comparisons are 
shown in table 1. Good agreement also can be expected at high supersonic and hypersonic 
crossflow Mach numbers where Newtonian theory by definition should be most applicable. 
The most doubtful regimes include the transonic crossflow Mach number regime and the 
supercritical Reynolds number regime. (These regimes are discussed further for circular 
cross sections in the following section.) 

Where reliable experimental crossflow drag data exist for a desired noncircular cross 

section, these data, of course, can be used. Then the values of Cj for the particular cross 

u n 

sections can be substituted in equations (2.20) and (2.21) in lieu of the product 

Cj ( C„/C„ ) , where Cj as now written is for a circular cross section only. Of course, 

a n n n o Newt » 

if the shape of the noncircular cross section varies along the body length, values of Cj must 

u n 

be substituted within the integral in the second term of equations (2.20) and (2.21), and a 
great deal of empirical input data may be necessary for some configurations. 

TABLE 1.- Cj AND CJC„ VALUES FOR TWO-DIMENSIONAL CYLINDERS OF 
u n o 

VARIOUS CROSS SECTIONS AT a = 90° AS COMPUTED BY NEWTONIAN 
THEORIES AND MEASURED AT SUBCRITICAL MACH AND REYNOLDS 
NUMBERS 


CROSS SECTION 

NEWTONIAN THEORY 
Cd n | Cn / Cn 0 

MOD. NEWT. THEORY 

™ w ,j 

Cd n | Cn / Cn 0 

MEASURED 

Udn | Cn / Cn 0 | REF. 

>o 

1 

1.33 1 1.00 
1 
i 

T 

1 

1.20 1 1.00 
1 
i 

1 1 

1.20 1 1.00 1 39 

1 1 

j i 

4 a/b- 2 

9 o/b-4 

i 

0.94 ] 0.50 
0.59 | 0.22 

L 

■■ \ 

0.85 0.50 

0.53 0.22 

i 

0.70 j 041 ] 40,41 
0.35 | 0.15 j 40 


r 

1 

1.65 1 1.75 

1 

i 

t 

1 

1.49 1 1.75 

1 
1 

1 I 

1 1 

1.60 1 1.89 1 41,42 
1 1 

1 1 

r-kw I* -0 - 0 

pp 1 k * 0.02 

tw k • 0.08 

Ll J k • 0.24 

k- 0.50 

2.00 i 1.33 
1.97 ] 1.33 
1.89 1.26 

1.68 j 1.14 
1.33 1 1.00 

i 

1.80 i 1.33 

1.78 | 1.33 

1.70 1.26 

1.51 | 1.14 

1.20 1 1.00 

i 

2.05 i 1.51 i 40 

2.00 | 1.48 | 41 

1.65 1.22 43 

1.12 | 0.85 | 43 

1.20 1 1.00 1 39 

i i 


NOTE: ALL C dn 's IN TABLE ARE BASED ON WIDTH OF 

CROSS SECTION, NOT EQUIVALENT d. 
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2.3 Empirical Input Values 


2.3.1 Crossflow drag coefficient- To compute C N and C m for the equations derived in 
the previous section, values of crossflow drag coefficient are needed for an “infinite 
length” or truly two-dimensional circular cylinder placed normal to an airstream. As pre- 
viously mentioned, is a function of both the Mach number and Reynolds number 
components normal to the cylinder longitudinal axis, and hence for a configuration at angle 
of attack it is a function of M n = M ^ sin a and R e n = Re sin a. In these simple relations, 
introduced previously as equations (2.3) and (2.4), M n is commonly called the crossflow 
Mach number and Re n , the crossflow Reynolds number. For circular cylinders, necessary 

“state-of-the-knowledge” plots have been prepared for the variation of Q with M and 

a n n 

Re n (figs. 1-3). 

Figure 1 shows the variation of C d with M n over the M fl range from 0 to 8. It was 
prepared from the data of references 44 through 49 and from data obtained recently by 
John M. Macha in the Ames 2- by 2-Foot Transonic Wind Tunnel. Also shown for reference 
are the theoretical variations predicted from Newtonian and modified Newtonian theories. 
Because of the close agreement of the Newtonian values with experiment at the higher Mach 
numbers, it is not surprising that computer programs utilizing Newtonian theories have been 
used successfully to predict space-shuttle-booster results in wind tunnels at hypersonic Mach 
numbers (see, e.g., ret. 50). Except for the transonic range, where data are very limited, the 
variation of Cj with M fJ is well documented in figure 1. 

In the transonic range, the black symbols in figure 1 represent values of C d obtained 

a n 

recently from pressure-distribution tests of circular cylinders of various diameters (1.9 to 
5.1 cm) at crossflow Reynolds numbers from about 1.3X1 0 s to 4.9X10 5 . Many values of 
C d n were initially computed from the extensive pressure distributions measured by Macha 
on the cylinders in the Ames 2- by 2-Foot Transonic Wind Tunnel. Because there was a 
general increase in C d with decrease in cylinder diameter d (but not Re n ) for M fJ values 
from about 0.9 to 1.2, plots ot C d vs. d were constructed, and values of Q were obtained 
by extrapolating the curves to d = 0. The black symbols in figure 1 represent these extrap- 
olated values of Cj , values that should come closest to representing data for no interfer- 
ence between the models and the wind tunnel. These data agree well with the rocket 
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flight-test results obtained in 1953 by Welsh (ref. 49), but these or similar-sized models 

should be tested further in a larger transonic wind tunnel. 

As shown in figure 1, there is a critical crossflow Reynolds number effect that can 

drastically lower the values of Cj at M n below about 0.5. For M less than about 0.5, if 

u n 

the crossflow Reynolds number Re n exceeds the critical value of about 2X10 S , the value of 

Cj decreases considerably. This variation is shown in greater detail in figures 2 and 3. 
a n 

Figure 2 gives the variation of C d with Re fJ for M n less than about 0.4. It has been 
well documented over the last 60 years (e.g., refs. 39, 43, 45, 51) that Cj^ = 1 .2 tor laminar 
boundary-layer flow and separation just before the critical Reynolds number of about 
Re n = 2X 10 s . At about Re n = 5X 10 5 there is evidence (e.g., refs. 52-54) of laminar 
boundary-layer flow around the front of the cylinder to an angular position of about 80° or 
90°, where the flow separates, undergoes transition, and reattaches at an angular position of 
about 110° to form a laminar separation bubble. Then the turbulent flow separates at some 
position downstream (an angular location of about 130 ). With a further increase in Rey- 
nolds number into the supercritical regime, the bubble decreases in size until the transition 
to turbulent flow moves upstream of the location of laminar separation, and the bubble 
disappears (ref. 54). From the low Cj value between about 0.15 and 0.30, C d increases 
gradually, at least for an increase in Re n up to about 5X10 6 . The supercritical Reynolds 
number regime has only been investigated recently in any detail (refs. 43 and 52-54), and 
there is still considerable uncertainty in the magnitude and trend of C d with Re n and M n . 
The shading in figure 2 indicates the approximate spread or uncertainty in C rf ^ based on 
known data. 

Jones, Cincotta, and Walker (ref. 54) probably have made the most detailed study of 

circular cylinders in supercritical flow. With the use of freon gas to obtain high they 

have shown that there is an effect of M n on the variation of C d with Re n . Figure 3 (taken 

from ref. 54) summarizes their Cj results for M n — 0.25 to 0.50. The reader is referred to 

u n 

reference 54 for their interpretation of these Cj results based on pressure-distribution and 

u n 

visual-flow studies. 

For noncircular bodies, as mentioned in section 2.2, experimental values of C d for 
noncircular instead of circular cross sections can be used with slight modification to equa- 
tions (2.20) and (2.21). These data are generally not available for crossflow Mach numbers 
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above critical. However, some data are available for subcritical crossflow Mach numbers, and 
these data should be used if the crossflow Reynolds number Re fJ exceeds the critical value 
(the Re n where C ^ drops rather drastically with a slight increase in Re n ). In table 2 some 

TABLE 2.- REFERENCES FROM WHICH EXPERIMENTAL VALUES OF C H VS 

a n 

Re n CAN BE OBTAINED FOR VARIOUS CROSS SECTIONS AND FLOW 
DIRECTIONS 
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references are listed from which experimental values of C ^ versus Re„ can be obtained for 

u n n 

various cross sections and flow directions. It should be noted that most experimental values 

of Q are based on cross-sectional width w and must be multiplied by w/d , where d is the 
u n 

diameter of the equivalent circular cross section. 

2.3.2 Crossflow drag proportionality factor— In the equations used to compute Cy 
and C m , r\ is the crossflow drag proportionality factor, that is, the ratio of the crossflow drag 
coefficient for a finite-length cylinder to that for an infinite-length cylinder. Cylinder drag 
coefficients from which values of 77 can be determined have been measured (to the author’s 
knowledge) only at very low subsonic Mach numbers (refs. 56 and 57). 

Values of 77 for circular cylinders at very low crossflow Mach numbers (from ref. 56) are 
plotted as a function of length/diameter ratio in figure 4. Values of 77 for flat plates are also 
plotted, but as a function of plate length/width ratio. The values for the flat plates are only 
slightly less than those for the circular cylinders. Thus, it is likely that figure 4 can be used 
to estimate values of 77 for many cross sections varying from circular to flat. However, these 
values may be acceptable only for very low crossflow Mach numbers. 

An indication of the variation of 77 with crossflow Mach number M n can be obtained 
by computing values of 77 from high-a C yy data (ref. 16) for slender bodies of revolution. 
From equation (2. 1 2), 

Cat ~ sin 2a cos (a/2) (A h /A r ) 

= — 7 T -/ Vw~S (2-22) 

n (AplA r )sm x a 

For two bodies of fineness ratio 10 and 12 (sketched in fig. 5), the variation of rjC j 

u n 

with M n (for M n = 0.4 to 1 .6) has been computed from equation (2.22) with the use of Cy 

data for values of a from about 45° to 60°. As shown in figure 5, the results for the two 

bodies agree closely. Now with the variation of Cj with M in figure 1 and the variation of 

a n n 

r]C with M n in figure 5, the variation of 77 with M n has been computed and the results 
plotted in figure 6. (The circular symbols denote the computed values over the M n range 
from 0.4 to 1.6; the square and diamond symbols represent values of 77 for very low 
crossflow Mach numbers obtained from fig. 4) 

For most supersonic and hypersonic values of M n , 77 probably can be assumed to be 
unity, an assumption indicated as being essentially correct from past investigations (e.g., 
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refs. 3, 6, 7, and 9). The greatest uncertainty appears to be within the transonic M n range, 
and here further research is desirable. In this study, however, figures 5 and 6 are used in lieu 
of better information. 


2.4 Formulas and Values of (C„/C„ ) and (C„/C„ ) for Various Cross Sections 

° SB l o Newt 

To use equations (2.20) and (2.21) for computing and C m , ratios of the section 

normal-force coefficients (C, 7 ) to those for the equivalent circular sections (C„ ) must be 

n o 

determined. The ratios ( C t JC„ ) from slender-body theory are used in the first term of 

o SB 

each equation, and the ratios (C„/C„ ) from Newtonian theory are used in the second 

o Newt 


term. 


Formulas of (C„/C„ ) 
o 

general cross sections encountered in missile and aircraft aerodynamics. Then, for several 


n ,~ n , and (C„/C„ ) are now presented for some of the more 
° SB ° Newt 


sample cross sections, values of (C„/C„ ) and (C M /C„ ) are plotted and compared. 

o SB n n o Newt 

2.4.1 Formulas of (C n /C„ ) - From slender-body theory (e.g., refs. 58-61), the 

n o SB 

ratio of C n for a winged-body cross section to that for the equivalent (same area) circular- 

body cross section can be determined for many cross-sectional shapes. In ref. 11, 

( C„/C n ) expressions are determined for winged-circular and winged-elliptic cross sections 
n n o SB 

(see sketches (c) ? (d), and (e)). 


S -*1 



(c) 




hH 



(e) 


Sketches (c), (d), and (e) 


For a winged-circular cross section with the wing planform perpendicular to the cross- 
flow velocity V n (sketch (c)), 



(2.23) 
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For a winged-elliptic cross section with the semimajor axis a and wing planform per- 
pendicular to the crossflow velocity V n (sketch (d)). 


where 



(2.24) 


k i = Oi 


(a + h) 2 
4aj 


a, ~ \ ( S + \A 2 + “ fl2 ) 


For a winged-elliptic cross section with the semiminor axis b and wing planform per- 
pendicular to the crossflow velocity V fl (sketch (e)), 


where 




+ b 2 ) 


k 2 


a 2 


(a + h) 2 
4a 2 



For an elliptic cross section without a wing (e.g., ref. 37), 


(2.25) 


( — | ^ cos 2 <p + — sin 2 0 (2.26) 

\C i 6 a 

\ . 

where 0 is the angle of roll about the body longitudinal axis, being 0° with the semimajor 
axis a perpendicular to the crossflow velocity and 90° with the semiminor axis b perpendic- 
ular to the crossflow velocity (see sketches (d) and (e)). 

2.4.2 Formulas of (CJCy. ) - From Newtonian impact theory, (C„/C„ 7 

n o Newt o Newt 

expressions also have been derived for winged-circular and winged-elliptic cross sections (see 
appendix A). 

For a winged-circular cross section with the wing planform perpendicular to the cross- 
flow velocity V fJ (sketch (c)), we obtain, from equation (A10), 
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© -Hh) 


(2.27) 


0 'Newt 

For a winged-elliptic cross section with the semimajor axis a and wing planform per- 
pendicular to the crossflow velocity V n (sketch (d )), we obtain, from equation (A 16), 


© 

' ° 'Newt 


3 

2 \ b 


- b 2 ja 2 


fl - (b 2 /a 2 )] 3 ' 2 


iog i (' + vn? ) 


I - (b 2 /a 2 ) a 


+ -- 1 


(2.28) 

For a winged-elliptic cross section with the semiminor axis b and wing planform 
perpendicular to the crossflow velocity V n (sketch (e)), we obtain, from equation (A21), 


© 


° 'Newt 


a 2 fb 2 


[( a 2 lb 2 )~ 1] 3/2 


tan* 1 




1 


{a 2 lb 2 ) - 1 b 


+ f -1 


(2.29) 


From Newtonian impact theory, an expression for (C_/C_ ) also has been derived 
f . . o Newt 

tor winged-square cross sections with rounded corners (sketch (0). From equation (B1 5) 


K s H /f»\ =l(i_ A) J r 


r = kw 


it)k 2 


0 < k < 0.5 (2.30) 


r 


L 

j 


where the equivalent diameter (from eq. (B12)) is 

|^1 - (4 - ir)k 2 


d = 2w 


(2.31) 


v n 


Sketch (f) 


2.4.3 Values of (C n /C n ) and (C n /C n ) - From equations (2.23) through 

° SB o Newt 

(2.29), values of (C n /C n ) and (C n /C H ) were computed for circular and elliptic 

° oB o Newt 

cross sections alone and with wings. The results are plotted and compared in figures 7 
through 10. 

In figure 7, the variation of ( C n \C n ) with axis ratio a/b is given for an elliptic cross 

section without wings. As previously noted in reference 9, values of (C /C ) from 

n ii/, 
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slender-body theory are reasonably close to those from Newtonian theory for many a/b 
values of interest. 

In figure 8, the variation of ( C fJ /C n ) with the ratio of wing semispans to body radius 

r is given for a winged-circular cross section. For s/r <2, the values from both theories are 

reasonably close, but with further increase in s/r, the values of (C n /C n ^) greatly exceed 

those of (CJC ) 

n n o Newt 

In figure 9, values of (C„/C„ ) are presented for a winged-elliptic cross section with the 

n ri Q 

semimajor axis a perpendicular to the crossflow velocity V n . For the axis ratios of a/b = 2 

and 3, the figure gives the variation of (C„/C„ ) with the ratio of semispans to semimajor 

' o 

axis a. As either a/b or s/a increases, the disagreement between the results from the theories 
increases. 

In Figure 10, values of (C„/C„ ) are presented for a winged-elliptic cross section with 

n n o 

the semiminor axis b perpendicular to the crossflow velocity V n . For axis ratios of a/b = 2 

and 3, the variation of (C„/C„ ) with s/b is given. There is closer agreement between the 

o 

values computed from the two theories for this cross-sectional arrangement than for the 

arrangement where the semimajor axis and wing are perpendicular to V n . 

Table 3, taken from reference 9, shows the variation of (C„/C„ ) with corner 

" n o Newt 

rounding k for square cross sections. The values of C n /C n computed from slender-body 
theory are reasonably close to those computed from Newtonian theory (eq. (2.30)). 


TABLE 3 - CJC„ FROM NEWTONIAN THEORY FOR SQUARE CROSS 

n n Q 

SECTIONS WITH ROUNDED CORNERS 


k 

0 

0.05 

0.1 

0.2 

0.3 

0.4 

0.5 

c„/c„ o 

1.33 

1.29 

1.25 

1.17 

1.11 

1.05 

1.00 


Corresponding values from slender-body theory vary from 1.19 at k = 0 (no corner 
radius) to 1.00 at k = 0.5 (completely circular cross section). 


2.5 Relative Influence of Crossflow Terms 

It is interesting to examine briefly the relative influence of the potential and viscous 
crossflow terms in the equations for Cyy and C m - For demonstration, Jorgensen (ref. 9) 
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compared the magnitudes of the terms for an ogive-cylinder body of fineness ratio 1 1 
(iflld = 5) at a = 0° to 1 80° and M ^ = 2.9. The computed values (by eqs. (2. 1 2) to (2.14)) 
are presented in figure 1 1 . For this body, the viscous crossflow term contributes most of the 
normal force at high values of a and, of course, all of the normal force at a = 90°. Although 
the slender-body potential term contributes relatively little to Cyy at high a, it has a signif- 
icant influence on C m . 

The results in Figure 1 1 are indicative of those computed for most slender bodies 
(ref. 9). However, as shown in chapter 5, the relative contributions of the crossflow terms 
can be modified considerably with the addition of thin lifting surfaces (wings and tails) to a 
body. 
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CHAPTER 3 

METHOD APPLIED TO BODIES OF REVOLUTION 

In recent years, there has been increased interest in the basic aerodynamics of slender 
bodies of revolution because of emphasis on achieving more maneuverability from missiles. 
Some designs that use thrust-vector control systems are being considered for missile flight at 
angles of attack ranging from 0° to 180° (e.g., refs. 13, 15, and 62). We will now assess the 
prediction method of chapter 2 by comparing predicted with measured longitudinal aero- 
dynamic coefficients tor various bodies of revolution. 

3.1 Cone-Cylinder and Ogive-Cylinder Bodies at = 2.9 

We will compare predicted longitudinal aerodynamic coefficients with those measured 
by Jemell (ref. 63) for a series of three cone-cylinder and four ogive-cylinder bodies at 
a = _ 5 ° to 180°. The free-stream Mach number is about 2.9 (Jemell quotes 2.86), and the 
Reynolds number based on body diameter is about 1 .25X10 5 . The bodies with various nose 
and aftersection fineness ratios were tested in the NASA-Langley Unitary Plan Wind Tunnel. 

Figure 1 2 shows the seven bodies considered along with values of the geometric param- 
eters required to compute the aerodynamic characteristics. For the cone-cylinder bodies 
(numbers 1-3), all of the geometric parameters are easily computed, but for the bodies with 
tangent ogive noses (numbers 4-7) the required values of A p /d 2 , V/d 3 ,x c /d, and A s /d 2 are 
not so easily obtained. Some convenient formulas (from ref. 9) for computing these param- 
eters for tangent ogives are given in appendix C. 

Equations (2.12) to (2.15), (2.18), and (2.19), along with the procedure outlined in 
chapter 2, have been used to compute the variation of C N , C A , C m and x ac /i with a for the 
seven bodies considered. The values of C a ^_ q0 and 180 ° used in e 9 uations (2-18) 

and (2.19) were computed in reference 9 with the assumption of turbulent boundary-layer 
flow. (Jemell (ref. 63) states that “boundary-layer transition was effected” by artificial trips 
throughout the a range.) The values include the contributions of fore pressure, base pres- 
sure, and the turbulent skin friction. The reader who wishes to make similar calculations is 
referred to reference 9. No attempt was made to estimate effects of wind-tunnel support 

interference. 
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In figures 13 to 16, computed values of Cyy, C A , and C m as a function of a are 
compared with the experimental results for the seven bodies. Generally, there is close 
agreement of the computed with the measured results, especially in the variation of Cyy and 
C m with a. As expected, because of the approximation formulas used to predict C A with a, 
the poorest agreement is between the predicted and measured values of C A . 

It is satisfying to find that effects of afterbody fineness ratio, nose fineness ratio, and 
nose shape on Cj y and C m are predicted so well. Generally the magnitudes of Cyy and C m 
increase with an increase in fineness ratio, just as the computed results predict. Figure 1 3 
shows the effect of afterbody fineness ratio for cone-cylinder bodies, all with fineness- 
ratio-3 conical noses. Likewise, figure 14 shows the effect of afterbody fineness ratio for 
ogive-cylinder bodies, all with fineness-ratio-5 ogival noses. Figure 15 shows the effect of 
nose fineness ratio for ogive-cylinder bodies, and figure 16 shows the minor effect of change 
in nose shape from conical to ogival for a given nose fineness ratio of 3. 

In figure 17, computed positions of aerodynamic normal-force center (symbols) are 
compared with measured positions (lines from ref. 63) for the seven bodies at M ^ = 2.9. 
The positions, x ac /t, are measured from the nose tip of each body in terms of the body 
length. As for the Cjy and C m results, the agreement of the computed with the measured 
values is reasonably close, especially for a near 90°. Note that symbols are used to denote 
computed values because only lines are given in reference 63 to denote the measured results. 

To assess the analytical method, plots of C m versus a may be omitted if plots of both 
C/V anc * x ac versus a are included. Hence, for conciseness, plots of C m versus a are omitted 
in the remainder of this study. Because precise prediction of C A versus a is beyond the 
scope of this study, plots of C A versus a are also omitted. 

3.2 Ogive-Cylinder Bodies at M — 0.6 to 2.0 

As shown in the previous section, the analytical method predicts the aerodynamic 
coefficients reasonably well for various cone-cylinder and ogive-cylinder bodies at a super- 
sonic Mach number of 2.9. We will now assess the method of predicting the C N and x QC 
results for four ogive-cylinder bodies at ot = 0° to 60° throughout the Mach number range 
from 0.6 to 2.0 (M^ = 0.6, 0.9, 1 .2, and 2.0). 
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Sketches of the four bodies considered are shown in figure 18. Bodies N 2 C X and N n C x 
were previously tested by Jorgensen and Nelson (ref. 16) in the NASA-Ames 6- by 6-Foot 
Wind Tunnel. Recently, these same bodies were retested in the same tunnel at = 0.6. 
They were also retested at = 0.6, 0.9, 1.2, and 2.0 with cylindrical extensions attached 
to the aftersections. Cylinder C 2 , two diameters long, was attached to body N 2 C X to form 
body N 2 C X C 2 , and cylinder C 3 , four diameters long, was attached to body (V 7 C, to form 
body N 7 C X C 3 . Thus, the bodies include nose fineness ratios of 2.5 and 3.5 (noses N 2 and 
N 2 ), aftersection fineness ratios of 7, 9, and 11 (aftersections C,, C X C 2 , and C,C 3 ), and 
overall fineness ratios of 9.5 (N 2 C X ), 10.5 ( N 2 C X ), 1 2.5 (N 2 C X C 2 ), and 13.5 (N 1 C X C 3 ). 

In figures 1 9-22, computed values of Cj y and (2 - x ac )/d as a function of a are com- 
pared with the measured results for the four bodies. The aerodynamic normal-force center, 
(2 - x ac )/d , is measured forward on each body from its base in terms of its body diameter. 
For bodies N 2 C X and /V 7 C, (figs. 19 and 20, respectively), the results for all Mach numbers 
are compared on the same plots, and one can observe, at a given value of c«, the variation of 
Cj y and ( SL-x ac )ld with M <*,. For the longer bodies, N 2 C X C 2 (fig. 21) and N 2 C X C 3 
(fig. 22), the results for each Mach number are compared on a separate plot. With these 
separate plots, there is less confusion in comparing predicted with measured values of 
(2 - x ac )/d. 

Generally, the variation of Cjy with a is predicted closely for each body throughout the 
r x and M 0 0 ranges considered. At the lowest Mach number, M m = 0.6, the (yy values previ- 
ously measured (ref. 16) for the shortest body, N n C x , are believed to be erroneous (fig. 20). 
The values from the retest agree closely with the predicted results. 

The aerodynamic normal-force centers are predicted best for the supersonic Mach 
numbers. Generally, for the subsonic Mach numbers, the predicted aerodynamic centers are 
more rearward on the bodies than the measured centers. For all Mach numbers and low 
values of a, the agreement of the predicted with the measured positions probably should be 
better than shown. The force and moment balance located inside each body was chosen to 
measure the large normal forces and pitching moments expected over the high a range. 
Thus, the accuracy of the experimental aerodynamic centers (determined from C m /C> y) at 
the low values of a is somewhat less than at the high values, and there is more scatter in the 
data. 
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In this study, we have concentrated primarily on predicting the static aerodynamic 
characteristics for bodies at high angles of attack. For bodies at low angles of attack, more 
theoretically exact prediction methods have been proposed. For example, for angles of 
attack less than about 20°, Perkins and Jorgensen (ref. 4) have shown that the agreement of 
theory with experiment at supersonic Mach numbers can be improved by replacing slender- 
body potential theory (e.g., the first term in eq. (2.10)) with Van Dyke’s hybrid theory 
(ref. 64) or Tsien’s linearized theory (ref. 65). They also suggest the use of an empirical 
modification to the viscous crossflow theory (e.g., the second term in eq. (2.10)). In this 
modification, an experimentally determined correlation curve for the distribution of the 
crossflow drag coefficient along the body length is used when there is laminar crossflow. It 
is questionable whether the extra computation needed in these more detailed methods is 
warranted for most engineering studies. However, if an accurate loading distribution over a 
body length is desired, a more detailed method should be considered (ref. 4). 

3.3 Predicted Effect of Change in Crossflow Reynolds Number 
from Subcritical to Supercritical 

At present there is a general lack of aerodynamic data for which predicted results can 
be compared for bodies of revolution at supercritical crossflow Reynolds numbers and 
subcritical crossflow Mach numbers. Jorgensen (ref. 9), however, has computed results 
which demonstrate what might be expected for a slender body of revolution. The body 
chosen for study consists of a tangent ogive nose of fineness ratio 5 with a cylindrical after- 
section of fineness ratio 6 (body 7 in fig. 12). 

For free-stream (and crossflow) Mach numbers less than about 0.4, the variation of Cm 
and x ac l% with a has been computed for free-stream Reynolds numbers of 10 s , 1 0 6 , and 
10 7 , and the curves are shown in figure 23. There is a significant effect of Reynolds number 
on both Cyy and * ac /£ throughout most of the a range. These curves, of course, reflect the 
strong influence of crossflow Reynolds number Re n on crossflow drag coefficient Cj for 

Tl 

two-dimensional circular cylinders (see fig. 2). 

As shown in figure 2, C j for a circular cylinder drops considerably as increases 

u n n 

from 10 s (subcritical) to 10 6 (supercritical), and then there is a gradual rise as Re n increases 
from 1 0 6 to 1 0 7 (in the supercritical range). There is much more uncertainty in the magnitude 
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of C,, at supercritical Re (such as 10 6 and 10 7 ) than at subcritical Re n (values less 

a n 

than about 2X10 5 ), and the shading in figure 2 indicates the uncertainty because of scatter 
in known data. 

In figure 23, the shaded bands in the C N and x ac /% curves for Re = 10 6 reflect the 
uncertainty in these curves resulting from the scatter in the data shown in figure 2. It is 
clearly evident, however, that this uncertainty in the curves is relatively small compared 
with the large effect of change in Reynolds number. 

Figure 23(b) shows the ratio of Cyy for the body at Re = 10 6 and 10 7 to for the 
body at the subcritical Re of 10 s . With this figure the effect of Re can be studied through- 
out the a range: for example, at a. near 90°, the body at Re = 10 6 develops only about 
25 percent of the C N developed at Re = 10 s , but at a < 10°, 100 percent of the Cpj is 
developed. Similar study of the movement of x flC /C with change in .Re can be made with the 
use of figures 23(c) and (d). 

Although no experimental data are available with which to compare these predicted 
effects for bodies of revolution, there are limited data for an early version of a noncircular 
space-shuttle body (ref. 66). Jorgensen (ref. 9) predicted some experimental trends (ref. 66) 
showing the decrease in C 'aj with a that results from an increase in crossflow Reynolds 
number from subcritical to supercritical. 
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CHAPTER 4 

METHOD APPLIED TO BODIES OF ELLIPTIC CROSS SECTION 

In this chapter, we compare predicted with measured normal-force and normal-force- 
center characteristics for three bodies of elliptic cross section at a = 0° to 60° and = 0.6 
to 2.0. First, however, we introduce the bodies considered, review the experimental test 
conditions, and present the equations used to compute Cy and C m for each body. 

4. 1 Bodies Studied and Tests at = 0.6 to 2.0 

Figure 24(a) shows the three bodies of elliptic cross section considered here. Planform 
views of the bodies as they were oriented (in five different configurations) for the tests are 
shown in figure 24(b). All of the bodies were tested previously and the results presented in a 
data report (ref. 1 8). 

The basic circular body B x consists of a circular-arc tangent-ogive nose of fineness 
ratio 3 followed by a cylindrical aftersection of fineness ratio 7. Bodies B 2 and B 3 have 
elliptic cross sections, and they have the same length and axial distribution of cross-sectional 
area as B x . Hence the fineness ratio of i/d = 10 for B x is also the equivalent fineness ratio 
for B 2 and B 3 , and all bodies have equal volumes. For B 2 , the ratio of the semimajor to the 
semiminor cross-sectional axis (a lb = 2) is held constant along the body length. Bodies B x 
and B 2 were investigated in 1 958 (ref. 6) only for a = 0° to 20° and M ^ = 2 to 4. Body B 3 
(new to the investigation in ref. 1 8) consists of the same nose shape as B 2 but has an 
afterbody section of variable a/b over four body diameters in length and a constant a/b = 2 
over the rear three body diameters (see fig. 24(a)). Photographs of B 3 in figure 25 enable 
one to establish a clearer mental image of this more complex body. 

Six-component aerodynamic force and moment coefficients were measured in the 
study of reference 1 8 for these bodies in the Ames 6- by 6-Foot Wind Tunnel. However, in 
this chapter, only the variation of Cjy and (C - jc ac )/d with a is considered. 

All bodies were tested at = 0.6, 0.9, 1 .2, 1 .5, and 2.0 and a = 0° to 58° . For the 
data used here, the Reynolds numbers, based on base diameter, are 6.5X10 5 at M <*, = 0.6 
and 0.9 and 3.8X1 0 s at M 00 = 1.2, 1.5, and 2.0. The elliptic bodies were tested at <p = 0° 
(flattest side of nose pitching against the flow) and 90° (see fig. 24(b)). 
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4.2 Equations Used to Compute Cjy and C m for Each Body 


For the general case of a slender body in which the cross-sectional shape varies along 
the length, the normal-force and pitching-moment coefficients can be computed from equa- 
tions (2.20) and (2.2 1 ): 
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The axial distance from the body base to the normal-force center (see sketch (g)) is 
then given in terms of the body diameter by 


(£ x ac ) _ C m + (£ x m^ 
d Cjy d 


(4.3) 



Sketch (g) 
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As previously discussed, (C„/C„ ) in equations (4.1) and (4.2) is the ratio of the 

n n Q SB 

local normal-force coefficient for the noncircular cross section to that for the equivalent 

(same area) circular cross section as determined from slender-body theory. The similar ratio 

( C t jC„ ) is determined from Newtonian impact theory. From equation (2.26), 

” n o Newt 



(4.4) 


where 0 is the angle of roll about the body longitudinal axis, being 0° with the semimajor 
axis a perpendicular to the crossflow velocity and 90° with the semiminor axis b perpen- 
dicular to the crossflow velocity (see sketches (h) and (i)). 



(h) (i) 

Sketches (h) and (i) 


From Newtonian theory, for the semimajor axis a perpendicular to the crossflow 


velocity V n , equation (2.28) reduces to 
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(4.5) 


For the semiminor axis b perpendicular to the crossflow velocity V n , equation (2.29) 
reduces to 
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For the bodies of this investigation, equations (4.1) and (4.2) can be simplified by 
integration. The simplified equations, which were used to compute C'yy and C m values for 
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each body, are presented next. In all the C ^ and C m equations, empirical input values of 17 

and Q from chapter 2 are used. 

u n 

4.2.1 Equations for bodies with constant a/b cross sections (bodies 5, and B 2 )~ Both 
bodies B x and B 2 have cross sections of constant a/b over the body length. For this 
condition, equations (4. 1 ) and (4.2) reduce to 
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4.2.2 Equations for body with variable a/b cross sections (body B 3 )- Body 3 has a 
midsection length of variable a/b, but the equivalent body of revolution is/?, for which 
(L4/dx = 0 rearward of the nose-cylinder junction. Thus, equations (4.1) and (4.2) reduce to 
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Some potential or inviscid normal force can be expected intuitively to be carried over 
past the nose shoulder of body 5 3 , even though the approximate slender-body term gives 
zero normal force over the aftersection. For body B 3 at 0 = 90° (fig. 24(b)), the span 
increases with length past the nose shoulder, and it might seem reasonable to compute the 
slender-body terms for Cyy and C m based on the maximum span (or base) geometry. For 
this estimation, 
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where 



4.3 Comparison of Computed With Measured Normal-Force and 
Normal-Force-Center Characteristics 

In figures 26 through 28, computed values of Cy and (£ - * ac )ld versus a are com- 
pared with measured values (from ref. 18) for the bodies studied. The comparisons are made 
for the bodies at a = 0° to 60° and free-stream Mach numbers of 0.6, 0.9, 1 .2, 1 .5, and 2.0. 

4.3.1 Bodies with constant a/b cross sections (bodies B x and B 2 )~ Generally, there is 
reasonably good agreement in figure 26 of the computed with the measured results for 
bodies B x and B 2 (bodies of constant a/b along the length). The agreement is, however, 
better at the supersonic Mach numbers than at the subsonic. These comparisons, along with 
previous successful comparisons in chapter 3, tend to validate the prediction method as a 
useful tool in body aerodynamic studies, at least for bodies with circular and elliptic cross 
sections of constant a/b. 

4.3.2 Body with variable a/b cross sections (body B 3 )— As for bodies B x and 5 2 , the 
predicted characteristics for body B 3 (with variable a/b) agree reasonably well with the 
measured results (fig. 27). However, the prediction of Cy generally is not as close as for bodies 
B x and B 2 . The prediction of (£ - x ac )/d at high a is remarkably close, especially at super- 
sonic Mach numbers. This close prediction of (8 - x ac )/d, however, might be somewhat 
fortuitous since the prediction of Cy is not nearly as close. 

The prediction of Cy with a is least accurate for body fi 3 oriented at 0 — 90°. Gen- 
erally, for this case, Cy is underpredicted at the lower values of a throughout the range. 
At the subsonic values of Cy is underpredicted over the entire a range studied. 

Note that these predictions were made with equations (4.9) and (4.10), which do not 
account for any potential or inviscid lift rearward of the nose section. Intuitively and 
analytically (ref. 64), some inviscid normal force can be expected to be generated rearward 
of the nose. This is especially true for body B 3 at 0 - 90° because the span increases with 
length past the nose, and it might seem reasonable to estimate the slender-body 
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potential-flow terms for C^y and C m based on the maximum spai? (see eqs. (4.11) and 
(4.12)). When the maximum span estimate is used, the prediction of Cyy witn a improves 
somewhat (see fig. 28). However, the prediction of (8 ~x ac )/d with a deteriorates, since the 
predicted values of (E - x ac )/d move well forward of the measured values (fig. 28). 
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CHAPTER 5 

METHOD APPLIED TO BODY-WING AND BODY-WING-TAIL CONFIGURATIONS 

In this chapter, we compare predicted with measured normal-force and normal-force- 
center characteristics for three body-wing configurations and three body-wing-tail configura- 
tions at a — 0° to 60° and = 0.6 to 2.0. First, however, we introduce the configuration 
models tested, review the experimental test conditions, and present the methodology used 
to compute Cyy and C m for each model. 

5.1 Configurations Studied and Tests at M <*, = 0.6 to 2.0 

Many body-wing and body-wing-tail model combinations have been tested, and the 
data are reported in references 20 and 21. Planform views of the model components tested 
in the many model combinations are shown in figure 29. These components include a basic 
circular body 5, , an elliptic body B 2 with a/b = 2, five flat-plate wings (IP, to W s ), and a 
combination horizontal and vertical tail T. The bodies B , and B 2 are two of the bodies con- 
sidered in chapter 4 (see fig. 24(a)). All the wings were designed to have the same planform 
area (16 d 2 ) if the wings extended into the body 5, to the axial centerline. Based on the 
phantom wing chord at the body centerline, the taper ratios for wings W X ,W 2 , and W 3 were 
0, 0.25, and 0.50, respectively (fig. 29(a)). They were also 0.25 for W 4 and W s (fig. 29(b)). 
Wings W t , W 2 , and W 3 (fig. 29(a)) had an aspect ratio of about 4; wings W 4 , W 2 , and W s 
(fig. 29(b)) has aspect ratios of about 5, 4, and 3, respectively. Pertinent planform dimen- 
sions of the exposed parts of the wings are given in tables in figure 29. The tail dimensions 
are given in figure 29(a). 

The tests (refs. 20 and 21) were conducted in the Ames 6- by 6-Foot Wind Tunnel at 
or = 0° to 58° and M ^ = 0.6 to 2.0. The Reynolds number, based on body diameter d , 
was about 4X10 5 . 

Results from these tests (refs. 20 and 21) showed that changing wing taper ratio from 0 
to 0.5 changed the aerodynamic characteristics very little. Also, changing wing aspect ratio 
from 3 to 5 changed the aerodynamic characteristics very little. Thus, to achieve more 
conciseness in the present investigation, only three body-wing configurations and three 
body-wing-tail configurations were considered. 
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Planform views of the configurations studied in the present investigation are shown in 
figure 30: the basic circular body B x with wings of aspect ratio 3 and 4 (fV 5 and W 2 , 
respectively) and the elliptic body B 2 with the wing of aspect ratio 3 (1P 5 ). These configura- 
tions (B x H/ s , B X W 2 , and B 2 W 5 ) were also investigated with the tail T attached (configur- 
ations B x W s T, B x W 2 T, and B 2 W 5 T). 


5.2 Methodology Used to Compute Cy and C m 
From equations (2.20) and (2.21), we obtain 
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The first terms represent the C N and C m values from slender-body potential-flow 
theory; the second terms represent the values from the viscous crossflow method modified 
by Newtonian theory. The first terms are not applicable, as written, for body-wing sections 
where the body cL4/dx values are zero or negative. Also, for body-wing and body-wing-tail 
configurations, more comprehensive methods from potential theory are available. 

For this study, the first-term (potential-flow) contributions to Cy and C m were com- 
puted from the linear method presented in NACA Report 1307 (ref. 38). This method, 
referred to as the N-K-P method (for Nielsen, Kaattari, and Pitts), is restricted to bodies of 
circular cross section with wings and tails that do not have swept-forward leading edges or 
swept-back trailing edges. It is further restricted to small angles of attack and small angles of 
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wing and tail incidence in which the forces are linear with angle. To obtain the wing-body 
interference, certain factors are defined that are the ratios of the lift on the components in 
combination to the lift on the wing alone. These ratios are obtained primarily by slender- 
body theory, but the wing lift is obtained by linear potential theory. Wing-tail interference 
is treated by assuming one completely rolled-up vortex per wing panel and evaluating the 
tail load by strip theory. 

To combine the N-K-P method with the crossflow method, the N-K-P potential terms 
must be multiplied by a correction factor (sin 2a)/2a to produce a more correct type of 
nonlinear behavior to these terms and to eliminate the potential contribution as a -*90°. 
With this modification, 
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Since the N-K-P method is restricted to bodies of circular cross section with wings and 
tails, a further assumption must be made to estimate potential theory values of Cyy and C m 
for noncircular bodies with wings and tails. The local widths of the noncircular body in 
planform are replaced by the local diameters of a circular body, thus keeping the overall 


wing and tail spans constant (see sketch (j))- The crossflow method, of course, requires no 


such assumption for noncircular bodies. 




Sketch (j) 



Assumed 
circular body 
for N-K-P 
method 
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For the way in which the crossflow method is formulated in chapter 2, values of 

crossflow drag coefficient C H for a circular cylinder are used. As shown in figure 1 , there is 

u n 

considerable variation of C., with crossflow Mach number M„ for values of M„ from about 

u fj fi rl 

0.4 to 3. Also, for M fl less than about 0.5, there is considerable change in Cj as the 
crossflow Reynolds number Re fl exceeds the critical value of about 2X1 0 s (figs. 2 and 3). 
These variations of Cj with M fl and R? n may be expected for near-circular bodies, but 

surely not for very flat bodies or winged bodies. From the data available, values of Q for 

a n 

flat bodies and plates do not appear to change nearly as appreciably with M n and Re n over 

the ranges shown in figures 1 and 2 (see, e.g., refs. 41, 66, and 67). Thus, for the body-wing 

configurations, it is likely that a constant value of C f i will give closer agreement of theory 

u n 

with experiment, especially at a where M fl is near or in the transonic regime. 

For flat-faced, two-dimensional configurations, reasonable values of crossflow drag 

coefficient can be computed from Newtonian or modified Newtonian theory (see, e.g., 

table 1 and ref. 66). For circular cylinders at low subsonic and hypersonic M , values of Cj 

' a n 

computed from Newtonian or modified Newtonian theory also agree reasonably well with 

experiment (see fig. 1). In this study, modified Newtonian theory is used to compute the 

circular-cylinder C , value that is substituted into equations (5.3) and (5.4). 
n 

From modified Newtonian theory, 

2 

C H = 7 - c n 

c n 3 Pstag 

= 1 .2 for C n =1.8 

Pstag 

ForM n > about 4, C p s(a ^ 18 trom perfect-gas relations (e.g., ref. 68). In this study for 

wing-body and wing-body-tail configurations, it is assumed that C H =1.2 for all values of 

a n 

M n (and hence, M oa ). 

The axial distance from the body base to the normal-force center is given (in body 
diameters) by 


® •* ac _ m W x m * 

d Cjy d 


(5.5) 


In this study, we present normal-force centers instead of C values. 
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The entire method has been computer programmed for Ames Research Center by 
Nielsen Engineering & Research, Inc., Mountain View, Calif. The program is similar to those 
given in reference 69, and it combines essential parts of the N-K-P and crossflow programs in 
reference 69. It is written in Fortran IV for the CDC 6600 or 7600 machines. No tapes, 
drums, or disks other than the standard input/output units are required. The running time 
for a typical case is less than 1 sec. 

5,3 Comparison of Computed With Measured Normal-Force and 
Normal-Force-Center Characteristics 

In figures 31 through 36, computed values of and (C - x ac )/d versus a are com- 
pared with measured values (refs. 20 and 21) for some of the body-wing and body-wing-tail 
configurations studied. As for the bodies alone, the normal-force coefficients are based on 
body cross-sectional area at the body base. Of course they can be easily converted to the 
more often used exposed wing planform area by dividing by A w /A^, where A w /A ^ = 1 5.92 
for the bodies with W s and 16.49 for the bodies with W 2 . The comparisons are made for the 
configurations at a = 0° to about 60°. For configuration B x W s , — 0.6, 0.9, 1 .5, and 2.0 

(fig. 31). For the other configurations (B x W 2 , B 2 W s , B x W s T x , B x W 2 T x , and B 2 W S T), 
M 0 o = 0.6 and 2.0. 

5.3.1 Body-wing configurations - Let us first look at the results in figure 31 for the 
circular body B { with the aspect ratio 3 wing W s . At M <*, = 0.6 (fig. 31(a)), the measured 
C j\j values are predicted closely by the modified N-K-P potential method only for a up to 
about 10° or 15°. Then this potential-flow method overpredicts Cyy for most of the higher 
a range considered. With even the potential-flow method overpredicting C N' the combina- 
tion of the potential-flow and crossflow methods greatly overpredicts Cyy. It is to be 
expected, however, that at a near 90° the combination method should predict C/y reason- 
ably well, since Cyy is given entirely by the crossflow method (modified Newtonian impact 
theory) at a = 90°. Despite the difficulty in predicting Cyy over most of the high a range, 
the positions of the normal-force center are predicted closely by the combination method. 

With increase in Mach number to ~ 0.9, the comparisons of computed with mea- 
sured results are not significantly improved (fig. 31(b)). In fact, the normal-force centers are 
not predicted as closely as at — 0.6. Not until the free-stream Mach numbers become 
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supersonic do the modified N-K-P and combination methods give reasonable estimates of 
Cyy. At Mao ~ 1-5, the measured Cy results are only underpredicted a small amount by the 
modified N-K-P method for a up to about 30° or 40° (fig. 31(c)). At Moo= 2.0, the Cy 
results are predicted best by the combination method throughout the a range studied, 
although the experimental results are still overpredicted by this method (fig. 31(d)). It 
appears that the combination method should improve in its ability to predict Cy as M^ 
increases throughout the supersonic range. 

With increase in wing aspect ratio from 3 to 4 (configuration B , W s to B x W 2 ), the 
comparisons generally are not changed significantly. As shown in figure 32(a), for = 0.6 
the modified N-K-P and combination methods still overpredict Cy significantly at high a. In 
fact, the overprediction for /?, W 2 is greater than for#, W s (compare figs. 31(a) and 32(a)). 
As for the lower-aspect-ratio configuration, the Cy y results are predicted best by the com- 
bination method as the Mach number is increased to — 2.0 (fig. 32(b)). 

When the body cross section is changed from circular to elliptical with the same aspect- 
ratio-3 wing (configuration B , W 5 to B 2 W 5 ), the Cy results for = 0.6 are predicted 
closely (possibly fortuitously) by the modified N-K-P method (fig. 33(a)). As before, when 
the Mach number is increased to M ^ = 2.0, the combination method still gives the best Cy 
prediction throughout the a range (fig. 33(b)). 

The break or significant deviation from linearity in the Cy curve with increasing angle 
of attack at the subsonic Mach numbers makes it extremely difficult to formulate a rational 
method for predicting Gy throughout the high a range. This break is attributed to flow 
separation over the wing upper surface. For wings of generally lower aspect ratio and higher 
leading edge sweep than those studied here, the break has been attributed to vortex break- 
down near the wing trailing edge of vortices shed from the leading edge. Vortex breakdown 
or bursting has been studied rather extensively by Wentz and Kohlman (ref. 70) for thin 
delta and modified delta wings with sweep angles from 45° to 85° at low speed. They have 
observed that as a increases the position of vortex bursting of the trailing vortices moves 
upstream toward the trailing edge and crosses the trailing edge at a specific a. Above this a, 
a loss of lift occurs on the wing due to vortex bursting, and the effect becomes progressively 
larger as a increases. Mendenhall and Nielsen (ref. 71) have more recently collected data 
from several investigators for the a value at which vortex bursting occurs at the trailing edge 
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of delta wings tested at low speeds. They were unable to correlate the data and suggested 
that the factors which control vortex bursting were not reproduced or controlled between 
the various sets of test data. The wings used in the present investigation were generally 
swept less than those studied by Wentz and Kohlman (ref. 70) and Mendenhall and Nielsen 
(ref. 71). However, a 45° delta wing similar to that of W x (fig. 29(a)) was investigated by 
Wentz and Kohlman (ref. 70). For this wing they failed to observe vortices, but they, of 
course, measured a loss in C ^ with increase in a over a particular value (near 20°). Despite 
the interesting research thus far, it seems that further research into the factors that control 
vortex bursting and flow separation is needed. The Cjy versus a data presented here indicate 
that this is especially desirable for subsonic Mach numbers. Some initial exploratory flow- 
field pictures from recent vapor-screen and oil-flow tests are presented in chapter 6. 

5.3.2 Body -wing-tail configurations- Now consider the results for the same body-wing 
combinations but with the tail added (figs. 34-36). Generally, the comparisons of computed 
with measured results indicate similar trends as for the body-wing configurations. At 
Moo = 0.6 and for the circular body B , , the C^y results over most of the investigated a range 
are underpredicted, even by the modified N-K-P method (figs. 34(a) and 35(a)). With change 
to the elliptic body, however, the Cy results (for B 2 W 5 T) are closely predicted up to a near 
50° (fig. 36(a)). At = 2.0, the Cy results are predicted best by the combination method 
(modified N-K-P plus crossflow), although as for the body-wing configurations the measured 
results are still overpredicted (fig. 36(b)). 

With the addition of the tail, there is generally more movement of the aerodynamic 
center (£ - x ac /d) with a at ~ 0.6 (compare, e.g., results for B x W s and B } W s T in 
figs. 31(a) and 34(a)). This movement, which takes place at a from about 10° to 40°, is only 
partially predicted by the combination method. It may be attributed to forebody and wing 
wake flow over the tail. Further investigation of the forebody and wing wake flow over the 
tail appears to be desirable. 

From the comparisons presented, it seems obvious that the methodology presented 
here represents only an initial step into the complex problem of predicting the aerodynamic 
characteristics of body-wing and body-wing-tail configurations to very high angles of attack. 
The reader interested in this field may wish to study several other initial approaches such as 
those of Mendenhall and Nielsen (ref. 71) and Axelson (ref. 72). One should also include the 
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Polhamus suction analogy for wings (ref. 73) and some of its various adaptations and exten- 
sions (e.g., refs. 71, 74, and 75). Much additional research is necessary in the high a field, 
and this research initially should include visual observations of the flow fields. In chapter 6, 
we will show some of the photographs obtained from an exploratory visual study of the 
flow over bodies alone, bodies with a wing, and bodies with a wing and a tail. 
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CHAPTER 6 

VISUAL OBSERVATION OF FLOWS OVER MODELS 

It is a well-established fact that as the angle of attack (or incidence) of a model is 
increased above about 5° or 10° some flow separates from the upper surface, and vortices 
and/or regions of How separation are formed. The basic vortex phenomena are easily 
described by referring to a body of revolution (see sketch (k)). At low to moderate angles of 
attack, two vortices are shed from the pointed body, and the vortex formation is sym- 
metrical (e.g., refs. 2, 5, 6, 13, 25, and 27). With further increase in angle of attack, the 
vortex formation may become somewhat asymmetrical. At some higher angle of attack, the 
feeding vortex sheets tear, and three or more vortices may appear (depending on the model 
geometry and free-stream flow conditions as well as angle of attack). Both symmetric and 
asymmetric vortex formations are illustrated in sketch (k). Although generally illuminating 


Symmetric vortex formation 




Asymmetric vortex formation and breakaway 
Sketch (k) 
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in concept, the illustrations in sketch (k) are somewhat simplistic, and further observations 
of vortex formations over many bodies alone and with wings and tails at various free-stream 
flow conditions are desirable. 

For observing the vortex formations over many models, the “vapor-screen technique” 
(e.g., refs. 2, 6, 76, and 77) is very useful and much less time-consuming than the more 
detailed pressure-probe-survey techniques (e.g., refs. 5 and 25). In this study, we have used 
the vapor-screen technique to observe the vortex formations over many of the same model 
configurations considered in chapters 3 through 5. In addition, we have added several ogival 
noses of different fineness ratio that attach to the body and body-wing-tail combinations 
previously described. To observe surface flow and especially flow separation positions, the 
oil-flow technique (e.g., refs. 6, 17, and 76) also has been used for many of the same 
configurations. 

In this chapter, photographs are presented that show the flows over the various bodies 
alone and in combination with a wing and a wing plus tail at a = 10° to 50°. The free- 
stream Mach numbers are 0.6, 0.9, and 2.0, and the Reynolds number is about 4.3X10 5 
based on body diameter. Before presenting the photographs, we will specify which models 
are considered and briefly review the vapor-screen and oil-flow techniques used for the tests 
in the Ames 6- by 6-Foot Wind Tunnel. 

6.1 Models Considered 

As in the previous chapters, the models considered are identified according to body B, 
nose N, cylinder C, cylinder strake S, wing W, and tail T. For the bodies of elliptic cross 
section (fig. 24), the roll orientation 4> is also specified. The configurations investigated and 
the figures in this report where the configuration dimensions are given are listed as follows: 


Configuration 

Figure 

fli = (V, c, 

24, 37 

N 2 c x 

37 

N 3 C t 

37 

a 4 c, 

37 

n 3 c,s 

37 
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Configuration 


Figure 


B 2 0 = 0° 

24 

B 2 0 = 90° 

24 

B 3 0 = 0° 

24 

B 3 0 = 90° 

24 

B,W S = NiCiWs 

29 

B x W 2 = N,C l W 2 

29 

B l W 2 T = N l C l W 2 

T 29 

n 3 c \ W 2 

29,37 

N 3 C t W 2 T 

29,37 


Plan-view sketches of these configurations are shown in figure 38. 

Figure 37, which has not been presented before, shows a series of circular-arc tangent 
ogive noses that attach to the circular cylinder C x . Noses N x , N 2 , and N 3 have fineness 
ratios of 3, 3.5, and 5, respectively. Nose yv 4 is formed by rounding the tip of a fineness- 
ratio-3.5 ogive (such as N 2 ) to give a resulting fineness ratio of 3. 

The thin wings and tail shown in figure 29 also attach to the cylinder C t , so many 
combinations of nose, cylinder, wing, and tail are possible. For this report only some 
representative combinations are considered. 

6.2 Vapor-Screen and Oil-Flow Techniques 

6.2.1 Vapor-screen technique - In the vapor-screen technique, the wind tunnel is run 
with moist air. In fact, water is added to the airstream as needed. 

For supersonic Mach numbers, as the moist air expands through the supersonic nozzle 
into the test section it cools, and the moisture condenses to form a fog. This fog is illumi- 
nated by a sheet of bright light produced by high-intensity mercury-vapor lamps and pro- 
jected through the tunnel window(s) and across the stream. This sheet of light appears as a 
uniformly lighted screen of fog particles in the absence of a disturbance. However, with a 
model in the stream the uniform distribution of fog particles is disturbed, and the model 
disturbance affects the light scattered by the water particles. Wakes and vortices typically 
appear as dark “holes” in the screen. 
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For very high subsonic — 0.9) and transonic Mach numbers, similar dark regions 
may appear at the vortex locations, but for lower subsonic flows (M^ = 0.6) light (conden- 
sation) regions often appear at the vortex locations. There are, in fact, cases where a light 
area will appear above a left wing panel and a dark area above the right panel, and both areas 
will indicate separated regions and/or vortex flow. The physics may not always be simply 
explained, but the reader interested in more detail is referred to the treatise of McGregor 
(ref. 77). 

Figure 39(a) shows a schematic drawing of the vapor-screen apparatus for the Ames 6- 
by 6-Foot Wind Tunnel. Two light-source boxes and various camera locations are indicated. 
In this study, two light boxes were used for observations at M = 0.6, but only one was 
required for necessary illumination at = 0.9 and 2.0. Each box contains six 900-W 
mercury-vapor lamps (BH-6) and mirrors to reflect the light through collimating slits. The 
light boxes (shown in fig. 39(b)) connect to a common shaft that passes over the tunnel, and 
the boxes can be pivoted in tandem and moved so that the light screen can cut the model 
longitudinal axis at the desired positions. A Honey well-Pentax spot light meter was focused 
on the vapor screen and moved with the light box shown in figure 39(b) to aid in maintain- 
ing consistent illumination for the photography. 

For the present investigation, a Hasselblad 70-mm still camera (model 500 EL/7) was 
mounted on the sting support rearward of the model base (fig. 39(c)). The camera was 
enclosed in a protective housing. The 16-mm gun camera shown mounted above the Hassel- 
blad in figure 39(c) was not used for this study. All models were painted black to minimize 
light reflection and to improve the quality of the photographs. The Hasselblad camera had a 
50-mm f4 (wide-angle) lens, and TRI-X (400 ASA) film was used. 

6.2.2 Oil-flow technique- Flow patterns on the model surfaces at angles of attack 
were visualized through use of the oil-flow technique. In this technique, the models were 
covered with a mixture of oil and titanium dioxide (Ti0 2 ) and then run wet in the wind 
tunnel. The formula for the mixture was 5 teaspoons of SAE 30 oil and 5 teaspoons of 
Ti0 2 , with about 3 drops of oleic acid added as an anticoagulant. To provide good contrast 
of the mixture with the models, all models were first painted a flat black. Photographs of 
the oil-flow patterns on the models were taken during each run with Hasselblad cameras 


46 



focused through a side window and a special small window at the top of the test section of 
the Ames 6- by 6-Foot Wind Tunnel. 


6.3 Photographs Obtained from Vapor-Screen and Oil-Flow Techniques 

6.3.1 Photographs from vapor-screen technique- Photographs taken with the vapor- 
screen technique for the models at = 0.6, 0.9, and 2.0 are shown in figures 40 to 54. 
Figure 40, presented as a prelude to the other figures, identifies some of the various items 
that appear in the other figures: the light sheet, light shadow, model support, model base, 
vortex regions, vortex feeding sheets, and local shocks. The two photographs in figure 40 are 
for body N 3 C l at a 30° and M ^ = 2.0. Light sheets are shown at the base of the nose 
(station 1) and at a more aft cylinder position (station 2). In the remaining vapor-screen 
photographs (figs. 41-54), a photograph of the flow field at the body base (station 3) is also 
included. The positions of the three flow-field stations are indicated on a model sketch in 
each figure. A model shadow similar to that shown in figure 40 usually does not appear on 
the photographs for M ^ — 0.6. As previously mentioned, at M ^ = 0.6, two light-source 
boxes were used in tandem, and light came from both sides of the wind tunnel. 

For a body of revolution (5, =A ^C, ) at M 00 = 0.6 (Fig. 41(a)), the simplistic flow 
model of the two rather tightly rolled-up vortices shown in sketch (k) is not evident, even at 
a < 30°. Rather there appear to be separation regions from both sides of the body that seem 
to almost coalesce along the lee side but do not roll up. In fact, two narrow separation 
sheets very close together appear to trail back from the nose over the body length. With an 
increase in a to about 40°, the separation regions from the nose appear more like the usual 
vortex regions. However, these vortex regions are very light in color in contrast to the dark 
regions at lower a, and these vortex regions become asymmetric as they trail back over the 
body. Also, more than two regions develop. 

With increase in Mach number to M = 0.9 and 2.0 (figs. 41(b) and (c)), the more 
traditional vortex formations appear, where there are two rather symmetric vortices (dark 
holes in the vapor sheet) shed from the body. Along both sides of the body, separation 
sheets “feed” the vortices, and, of course, the vortices grow with movement from the nose 
to the base of the body. At M ^ = 2.0, local shocks from the vortex regions appear 
(fig. 41(c)). 
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From observation of figures 41 to 43, the effect of increasing nose fineness ratio 
from 3 (nose N t ) to 5 (nose ;V 3 ) can be studied. Generally, with an increase in nose fineness 
ratio, the vortex formation becomes asymmetric at a lower a. Compare, for example, 
photographs forA/,C, at a % 30° (fig. 41 ) with those for N 3 C, at a as 30° (fig. 43). For all 
of these bodies, as for the remaining models to be discussed, vortex asymmetry is the worst 
at the lowest Mach number investigated = 0.6). Asymmetric flow separation and 
vortex asymmetry are accompanied by undesirable side forces (considered in chapter 7). 

When the tip of a fineness-ratio-3.5 nose (,V 2 ) was rounded to make a fineness-ratio- 
3 nose (N 4 ), a strange vortex pattern developed at M = 0.6 for the body A 4 C, at a ^ 40° 
and 50 . This pattern (fig. 44(a)) consists ot two very symmetric vortices from the nose 
located above two separation regions stacked one on top of the other (at station 2). The 
entire unusual pattern appears to be symmetric even up to o; ~ 50°. However, at the higher 
Mach numbers (A/^ = 0.9 and 2.0), the patterns are again similar to those for A^Cj (com- 
pare, e.g., figs. 42 and 44). 

When strakes were attached to the side of the cylinder C, of body W 3 C, , the vortex 
patterns became more symmetrical for a up to about 30° (compare figs. 43 and 45). How- 
ever, at about 40° and above, there was no effect. Apparently, the asymmetric pattern, 
originating with the fineness-ratio-5 nose, could not be influenced by the strakes back on 
the cylinder. Note that the hand of the asymmetry of the vortices can be either left or right, 
and occasionally the pattern will switch even while it is being observed during a test run. 

When the body cross section was changed from circular (body 5, ) to elliptic (body B 2 
with ajb ~ 2), the vortex patterns became more symmetric with the body oriented at 0 = 0° 
(flattest side toward the flow). This can be seen by comparison of figures 41 and 46. How- 
ever, when the elliptic body B 2 was rolled to 0 = 90°, the vortex patterns became more 
asymmetric (see fig. 47). 

From tests ot body B 3 , the body of elliptic cross section with variable ajb, the vortex 
symmetry was influenced mostly by the nose. As shown in figure 48, when B 3 was oriented 
at 0 = 0° (flattest side of nose toward the flow), the vortex patterns were essentially 
symmetric at all test conditions. However, when B 3 was rolled to 0 = 90° (thinnest side of 
nose toward the flow), the patterns became quite asymmetric at many test conditions (see 
fig. 49). 
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When wing W 5 of aspect ratio 3 was attached to body Z?, , the body vortex growth over 
the wings was retarded, but extensive regions of flow separation above the wing appeared 
(fig. 50). These regions were especially large and diffuse for - 0.6 and 0.9 (figs. 50(a) 
and (b)). For = 2.0, distinct vortices formed, originating at the forward wing-body 
juncture at the leading edge, on the wing upper surface, and at the wing tips. Thus, at the 
base of the body there were a pair of vortices from the nose, a pair from the wing-body 
juncture, a pair from the wing leading edge, and a pair from the wing tips (see fig. 50(c)). 

In contrast with the flow at M = 2.0, the flow at M = 0.6 is very diffuse, and there is 
the possibility that the phenomenon of vortex bursting or breakdown (e.g., refs. 70 and 7 1 ) has 
taken place at the higher angles of attack. This phenomenon (mentioned in chapter 5) might 
be indicated in some of the vapor-screen pictures by the light-colored diffuse vortex regions. 
Such regions can be observed at = 0.6 for both bodies alone and with wings, especially at 
the higher angles of attack. The phenomenon of vortex bursting might help explain the fact 
(demonstrated in previous chapters) that the Cyy characteristics for the bodies alone and with 
wings and tail are generally predicted best at the supersonic Mach numbers. The prediction 
methods, of course, do not account for this phenomenon. 

When wing W 2 of aspect ratio 4 was attached to body fij , similar flow patterns were 
obtained as with wing W s of aspect ratio 3. This can be confirmed by comparing the 
photographs in figures 50 and 51. Note, however, that the vortices shed from the wing tips 
offl, W 2 at M oq = 2.0 lie outside the photograph frames. 

With the addition of the tail T to configuration B l W 2 , there were no appreciable 
changes in the flow patterns (see photographs in figs. 51 and 52). 

When wing W 2 was attached to body A r 3 C|, the nose-cylinder configuration that 
produced the greatest asymmetry at M ^ = 0.6 and 0.9, the vortex asymmetry from the 
fineness-ratio-5 ogival nose still persisted (compare figs. 43 and 53). This asymmetry also 
persisted when tail T was attached to the configuration (see fig. 54). 

6.3.2 Photographs from oil-flow technique — To obtain visual indications of the flow 
over the surtace of the models, the oil-flow technique was used. Photographs from this 
technique support those from the vapor-screen technique in that flow-separation positions 
arc quite clearly defined. As illustrated in a schematic of a body crossflow plane (sketch (1)), 
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oil ridges form near separation regions, and How symmetry or asymmetry can be correlated 
with results from the vapor-screen technique. 



Sketch (I) 

Photographs taken with the oil-flow technique for some selected models at M ^ = 0.6, 
0.9, and 2.0 are shown in figures 55 to 65. Both planform and side views are shown at the 
specified angles of attack. To avoid a tunnel support at the top of the test section, the 
planform views were taken with the camera somewhat off center, so the model planforms 
are not completely symmetrical. For convenience in comparing the oil-flow and vapor- 
screen photographs, the figures for each technique are indexed as follows: 


Configuration 

Oil-flow 

figures 

Vapor-screen 

figures 

B i = /V,C, 

55 

41 

N i C l 

56 

42 

n 3 c, 

57 

43 

yv 4 c, 

58 

44 

n 3 c\s 

59 

45 

B 2 0 = 0° 

60 

46 

B 2 0 =90° 

61 

47 
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Configuration 

Oil-flow 

figures 

Vapor-screen 

figures 

o 

O 

II 

05 

62 

48 

B 3 4> = 90° 

63 

49 

B { W 2 T = N l C l W 2 T 

64 

52 

n 3 c\ w 2 t 

65 

54 


The reader, of course, can make detailed observations of the comparable oil-flow and 
vapor-screen photographs for configurations of particular interest. Generally, when there is 
asymmetric flow separation from a model surface (indicated by the oil flow), there is also 
asymmetric arrangement of the vortices in the flow field (indicated by the vapor screen). 
Again, note that the hand of the asymmetry can be either left or right. Thus, because the 
oil-flow and vapor-screen photographs were not taken simultaneously, it is possible to 
observe a left-hand separation asymmetry and a right-hand vortex asymmetry or vice versa. 
The observer also should be aware that some local flow disturbances in the oil flow result 
from joints and wax-filled screw holes in the bodies. 
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CHAPTER 7 

EXPERIMENTAL SIDE FORCES ON MODELS AT 0 = 0° 

When models are pitched to high angles of attack, side forces can occur on the models 
even at zero sideslip angle (e.g., refs. 1 2-21 , 78, and 79). These side forces generally occur at 
angles of attack between about 20° to 60° and in the subsonic-transonic Mach number 
range. They result from asymmetric flow separation and vortex flow over the leeward side 
of the models as shown in the photographs presented in chapter 6. 

As noted in a recent paper by Keener, Chapman, and Kruse (ref. 78), some aircraft 
have been lost due to uncontrolled flight at high angles of attack, and some of the loss in 
controllability might have originated from undesirable side forces and yawing moments 
attributed to flow separation and vortex asymmetry. Research on this phenomenon recently 
has increased considerably because the flight envelopes of modern aircraft and missiles are 
being extended into the higher angle-of-attack range. 

We will now present and discuss some side-force data obtained recently (refs. 16, 
18-21) in the Ames 6- by 6-Foot Wind Tunnel for most of the model configurations 
considered in previous chapters (especially chapter 6) and shown in figure 38. For the bodies 
alone, we will discuss the effects of nose-fineness ratio, Mach number, nose-tip rounding, 
afterbody side strakes, and elliptic cross section. Then we will consider the effect of adding 
a wing and a wing plus tail to a body. 


7.1 Bodies Alone 

7.1.1 Effects of nose fineness ratio and Mach number - In figure 66, the effect of 
nose fineness ratio Cyy/J on side-force coefficient Cy and center position (C - x s f)/d for 
ogive-cylinder bodies of revolution is shown. Both Cy and (C - x^jd, measured from the 
body base, are plotted as a function of angle of attack a for a up to 60°. Plots are presented 
for A/qq = 0.6, 0.9, and 1.2. The models were tested at Mach numbers up to M ^ = 2, but 
there were no side forces above about 1.2. The magnitudes of the side forces decrease 

with increasing Mach number. 

It is readily apparent from figure 66 that the largest side forces were obtained with the 
fineness-ratio-5 nose 7V 3 attached to the cylinder afterbody C x (of fineness ratio 7). With 
decrease in nose fineness ratio, the side forces decrease, and they almost disappear over most of 
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the a range for the body with the fineness-ratio-3 nose. For body N 3 C i (with the fineness- 
ratio-5 nose) at = 0.6, it has been shown that the maximum value of side force can 
become about 40 percent of the normal force (ref. 16). Keener et al. (ref. 78) have recently 
shown that the maximum value of Cy increases even more (to the order of Cy or greater) as 
the Mach number is decreased from 0.6 to 0.25. There is no question that the side-force 
coefficients increase with decrease in Mach number and/or increase in nose fineness ratio. 
Wardlaw and Morrison (ref. 79) also support these conclusions from their recent correlations 
of collected data. 

Note that the side-force centers seem tc start well back on the cylinder; then as a 
increases above about 25°, they move forward onto the nose (see, e.g., fig. 66(a)). At some 
much higher a (say about 45° to 55°) they then tend to move back onto the cylinder. Note 
also that for A l 3 C } , the body with the largest side forces, the maximum values of Cy are 
located well forward on the body. 

The side forces can be studied in conjunction with the vapor-screen photographs shown 
in chapter 6. In figure 67, the data for N 3 C 1 at = 0.6 and 0.9 are plotted along with the 
vapor-screen photographs (from figs. 43(a) and (b)) taken at a = 10°, 20°, 30°, 40°, and 
50°. The vapor-screen photographs shown in figure 67 were taken at a crossplane station 
3.5 diameters forward of the model base. It is obvious that when the vortices become 
asymmetric (between a = 20° and 30°) the side forces develop. So long as the vortices 
remain steadily asymmetric, the side forces remain. Of course, as shown in chapter 6, an 
asymmetric vortex system results from asymmetric flow separation from the body surface. 
There is obviously a need to analytically model the instability process that leads to asym- 
metric separation, asymmetric vortices, and undesirable side forces. 

In lieu of a theoretical method for computing the angle of attack for onset of side 
force, experimental procedures must be used. Keener- et al. (ref. 78) have found that the 
onset angles can be roughly correlated with nose fineness ratio (or semiapex angle) as 
families of curves of constant afterbody. They have found that the angle of onset of side 
force for a given body is essentially invariant with Mach number, and, for a nose with no 
afterbody, the onset angle is given approximately by two times the semiapex angle. The 
general effect of increasing afterbody length is to decrease the angle of onset; that is, the 
longer the afterbody the smaller the angle at which a side force is first encountered. A 
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somewhat similar finding is reported by Wardlaw and Morrison (ref. 79). They also have 
produced (ref. 79) a rough correlation of data for the angle of attack at which the maximum 
side force is observed. It tends to decrease with increasing Mach number and body fineness 
ratio. 

Keener et al. (ref. 78) also have made a rough correlation of data for the “upper-limit” 
angle of attack at which the static side force disappears and the wake flow becomes essen- 
tially oscillatory like a Karman vortex street from a two-dimensional cylinder. This upper- 
limit a also tends to decrease with increasing Mach number, varying from about a maximum 
of 80° at M x = 0.25 to a minimum of 50° at M = 0.9. 

7.1.2 Effect of nose-tip rounding — When the tip of a fineness-ratio-3.5 nose (N 2 ) was 
rounded to make a fineness-ratio-3 nose ((V 4 ), the side forces at M ^ = 0.6 and 0.9 almost 
disappeared. This is shown in figure 68 where Cy and its center position are plotted against 
a for bodies N, C, , N 2 C , , and /V 4 Cj . As discussed in chapter 6, this nose rounding brought 
more symmetry to the flow field. The round-nosed body JV 4 C , , however, appears to be no 
better than the sharp-nosed body N, C x of the same fineness ratio ( Hy/d = 3). 

7.1.3 Effect of afterbody side strokes When strakes were attached to the sides of the 
cylinder C x of body N 3 C X (K/y/d = 5), the side forces were not significantly changed. As 
shown in figure 69, the variation of Cy with a was changed somewhat, but the maximum 
values of Cy were about the same. 

7.1.4 Effect of elliptic cross section - The effect of elliptic cross section and roll angle 
on side-force coefficient and position is shown in figure 70. Results are compared for body 
B x (a circular body of i/d= 10) and body B 2 (the equivalent elliptic body of constant 
alb — 2) at = 0.6 and 0.9. With the elliptic body B 2 oriented at 0 = 0° (flat side 
pitching against the free-stream flow), the side-force coefficients are very small and close to 
those for B x . However, when B 2 is rolled to 0 = 90°, the side-force coefficients increase 
considerably. According to reference 18, Cy becomes more than twice C y at a = 50° for 
B 2 at 0= 90° and = 0.6. As shown in chapter 6, the separation and vortex patterns were 
very symmetric for B 2 at 0 = 0°, but they became very asymmetric for B 2 at 0 = 90°. 

From tests ol body B 3 , the body of elliptic cross section with variable alb , it was 
found that the wake flow-tield asymmetry and side forces were influenced mostly by the 
nose. As shown in figure 71. when B 3 was oriented at 0 = 0° (flattest side of nose pitching 
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against the free-stream flow), the side-force coefficients were very small and close to those 
for B x (the equivalent body of revolution). However, when B 3 was rolled to <p = 90°, 
undesirable side forces developed, and, as shown in chapter 6, the vortex flow was quite 
asymmetric (fig. 49). 

As for the bodies of revolution, the maximum values of side-force coefficient decreased 
as the free-stream Mach number increased. At M 00 = 1.5 there were essentially no side 
forces measured on any of the elliptic bodies, even for the worst roll orientation (ref. 1 8). 

7.2 Body-Wing and Body-WingiTail Configurations 

7.2.1 Effects of adding a wing and a wing plus tail to a body - In figure 72, the effects 
on side-force coefficient and side-force position of adding a wing and a wing plus tail to a 
body are shown. The body N 3 C\, which has a tineness-ratio-5 nose, produced the largest 
side forces for the bodies ot revolution tested. As shown in figure 72, these side forces and 
their positions remain about the same with the thin wing of aspect ratio 4 attached or 
even with the tail T added. (Dimensiohs for the wing and tail are given in fig. 29 and for the 
body in fig. 37.) 

From these comparisons it can be concluded that the most important influence comes 
from the body nose. As mentioned in chapter 6, the vortex asymmetry that appeared from 
the fineness-ratio-5 nose, when tested with only the afterbody cylinder C, , still persisted 
when the wing and the wing plus tail were attached (see e.g., figs. 43, 53, and 54). 

7.2.2 Effects of wing aspect ratio and taper ratio- For thin wings of about equal 
planform area (fig. 29) but with aspect ratios of 3, 4, and 5, there was essentially no effect 
of aspect ratio on the measured side forces for the body-wing models tested (ref. 20). 
Likewise, a change in taper ratio from 0 to about 0.5 resulted in no appreciable side-force 
effect (ref. 20). The results were also unchanged when the tail T (fig. 29) was attached to 
the wing-body models (ref. 21). 
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CHAPTER 8 

CONCLUDING REMARKS 


A review and an extension of an engineering-type method have been presented for 
computing the normal-force and pitching-moment coefficients for slender bodies of circular 
and noncircular cross section alone and with lifting surfaces. In this method, a semiempirical 
term representing viscous separation crossflow is added to a term representing potential- 
theory crossflow. Jn the generalized equations written for C'yy and C m , ratios are required of 
the local normal-force coefficient per unit length for the cross section of interest to that for 
the equivalent (same area) circular cross section. These ratios are given both from slender- 
body and Newtonian theories. Formulas and numerical values of these ratios are included 
here for winged-elliptic and winged-square cross sections, the square cross sections having 
rounded corners if desired. 

In computing normal-force and pitching-moment coefficients for the bodies alone, 
slender-body theory was used for the term representing potential crossflow. In computa- 
tions for the bodies with thin wings and a tail, the linearized potential-flow method of 
Nielsen, Kaattari, and Pitts was used, modified for high angles of attack. 

For many bodies of revolution, computed aerodynamic characteristics were found to 
agree well with measured results for investigated free-stream Mach numbers from 0.6 to 2.9. 
The angles of attack ranged from about 0° to 180° for = 2.9 and from about 0° to 60° 
for = 0.6 to 2.0. Agreement of predicted with measured results was best at supersonic 
Mach numbers. 

For several bodies of elliptic cross section, measured results were also predicted reason- 
ably well over the investigated Mach number range from 0.6 to 2.0 and at a = 0° to 60°. As 
for the bodies of revolution, the predictions were best for supersonic Mach numbers. The 
predictions were better for a body in which the cross-sectional shape (alb) remained con- 
stant over the length than for a body in which it varied. Although the prediction technique 
probably can be improved with further research, it is felt that the predictions are suffi- 
ciently accurate for most preliminary design studies. 

For body-wing and body-wing-tail configurations with wings of aspect ratio 3 and 4, 
measured normal-force coefficients and normal-force centers were predicted reasonably well 
at the upper test Mach number of 2.0. However, with a decrease in Mach number to 
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Moo = 0-6, the agreement for Cjy rapidly deteriorated, although the normal-force centers 
remained in good agreement. At the subsonic Mach numbers (M^ = 0.6 and 0.9) and angles 
of attack above 10° or 20°, the measured results were even overpredicted somewhat by just 
the potential-flow term of the combination method. 

From vapor-screen and oil-flow studies, it was observed that the flow was completely 
separated over the wings at the subsonic Mach numbers for a > about 20°. There was some 
evidence of vortex bursting or breakdown for some of the models (including bodies alone) 
at subsonic Mach numbers. At — 2, however, discrete vortices from the bodies and 
wings were observed at the base of the models. 

For many body, body-wing, and body-wing-tail configurations, vapor-screen and oil- 
flow photographs were obtained for = 0.6, 0.9, and 2.0 and a = 10°, 20°, 30°, 40°, and 
50°. It has been observed that, when the separation and vortex patterns were asymmetric, 
undesirable side forces could be measured on the models even at zero sideslip angle. These 
side forces generally originated when the angle of attack exceeded about 20°. 

For bodies alone, the side forces can be significantly affected by changes in Mach 
number, nose-fineness ratio, nose bluntness, and elliptic cross section. The side-force coeffi- 
cients decrease with increase in Mach number through the subsonic-transonic range and 
disappear with increase in Mach number into the supersonic range. 

From tests (at 0.6 < < 2.0) of tangent-ogive noses connected to a circular cylinder 

of fineness ratio 7, it has been found that the side-force coefficients increase from about 
zero for a fineness-ratio-3 nose to a maximum of about 40 percent of the normal-force 
coefficient for a fineness-ratio-5 nose. Other researchers have reported even more increases 
in Cy (of the order of Cy or greater) as the Mach number is decreased from 0.6 to 0.25. 

Nose-tip rounding significantly decreased the side-force coefficients for a circular body 
with a sharp-nosed ogive of fineness ratio 3.5. However, the beneficial decrease was no 
greater than that obtained by merely using a sharp-nosed ogive nose of the same fineness 
ratio (fineness ratio 3) as that for the resulting blunted nose. 

When strakes were attached to the sides of the afterbody cylinder of an ogive-cylinder 
model with an undesirable fineness-ratio-5 nose, the side forces were not significantly 
changed. Thus, the influence of the nose was dramatized. 
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Changes in cross section from circular to elliptic for an ogive-cylinder body (of 
fineness-ratio-3 nose and fineness-ratio-7 afterbody) produced some interesting effects on 
side-force coefficient. At all test Mach numbers (0.6 to 2.0) and angles of attack (0° to 60°), 
the side-force coefficients were generally small or negligible for the circular body and the 
generated elliptic bodies at 0 = 0° (flat side pitching against the stream crossflow). With the 
elliptic bodies at 0 = 90°, however, some values of side-force coefficient became as large as 
twice the values of normal-force coefficient at the same high angles of attack. From the 
standpoint of reducing undesirable side forces at high angles of attack, it was always found 
best to have the flattest side of the elliptic-body nose pitching against the stream crossflow, 
even when alb was not constant over the body length. 

Undesirable side forces measured for body-wing and body-wing-tail models were gen- 
erally about the same as those measured for the bodies alone. As for the bodies alone, the side 
forces developed at subsonic Mach numbers for a > about 20°. Also, as for the bodies alone, 
the side forces and yawing moments increased with increase in nose fineness ratio. Fineness 
ratios greater than 3 produced the largest side forces. No effects of wing aspect ratio or 
taper ratio were observed. From these comparisons, it can be concluded that the most 
important influence comes from the body nose. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, California 94035, March 18, 1977 
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APPENDIX A 


DERIVATION OF (C„/C„ ) FOR WINGED-CIRCULAR 
" n o Newt 

AND WINGED-ELLIPTIC CROSS SECTIONS 


Winged-Circular Cross Section with Wing Planform Perpendicular 
to Crossflow Velocity 


For blunt configurations of the type 
shown in sketch (m), it is assumed that each 
elemental particle of fluid strikes the 
configuration at velocity V n and thereupon 
loses its normal component of momentum. 
This leads to the well-known Newtonian 
expression for the pressure coefficient: 



Cp = 2 sin 2 5 (Al) 



Sketch (m) 


where S is the local angle that a tangent to a forward-facing surface makes with the 
free-stream direction (sketch(m)). In Newtonian flow, the pressure coefficients over the 
rearward face are assumed to be zero. The total section crossflow drag coefficient (based on 
body diameter, d = 2 r) is then given by 


c d = 



(A2) 


For a circular section, as shown in sketch (m). 


x 2 + y 2 = r 2 


(A3) 

c. dy x 

tan o = -r- = - — = — 

<±v y 


(A4) 
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and 


sin 2 5 -- 


tan 2 5 


= 1 -=^ 


.,2 


1 + tan 2 5 

Equations (A1 ) and (A5) are substituted into equation (A2) to give 


C d = ~ 


7 f 2 (' - $) iy + 7 f liy 


since 


C n = 2 sin 2 6=2 
rwing 


From integration of equation (A6), we obtain 

Q = 7 +2 (7- >)= 2 (y 

where, for the equivalent circular section, 


c = i 

3 


(A5) 


(A6) 


(A7) 


(A8) 


(A9) 


Thus, the ratio of the normal-force coefficient per unit length for the winged-circular 
section to that for the equivalent circular section is given by 


/ C n \ _ ( C<1 \ = 3 A 

\ C n / \ C d j 2 b 3/ 

' ° Newt ' n ^Newt 


(A10) 


Winged-Elliptic Cross Section with Semimajor Axis and Wing Planform 
Perpendicular to Crossflow Velocity 

For a winged-elliptic cross section where the semimajor axis a and the wing planform 

are perpendicular to the crossflow velocity V n (see sketch (n)), the same procedure is 

followed in the derivation of (C^/C„ ) as introduced for the winged-circular cross 

° Newt 

section. Both the basic equations, (Al) and (A2), are used again. 
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For an elliptic body, as shown in sketch (n). 



^ + iL 2 =l 

b 2 a 2 


(All) 




i 


and 

sin 2 5 = 


2 r _ tan 2 5 


a 2 - y 2 


1 + tan 2 5 a 2 - y 2 [ 1 - (b 2 /a 2 )] 


(A 13) 


Sketch (n) 


With equation (A13) combined with equations (Al), (A2), and (A7), we obtain the 
section crossflow drag coefficient: 


c = 4 f a a 2 - >> 2 
d d J a 2 - y 2 [ i - (b 2 /a 2 )] 


d y + 


r 


dy 


(A14) 


o ~a 

where d is the diameter of the equivalent circular body and is given by 

d = 2 sjab 


(Al 5) 


Also, as for the winged-circular cross section, the crossflow drag coefficient for the equiva- 
lent circular cross section is given by Cj = 4/3. 

a n 

Thus, by integrating equation (A 14), substituting 2 sjab for d, and dividing the result 
by 4/3, we obtain 


6) 

' °\Vewt 


±y* 

2 V b 


- b 2 ja 2 


[1 - (b 2 /a 2 )] 3 ' 2 


log 


[f (' + f 7 ^) 


4- 


! + L 

1 - (b 2 /a 2 ) a 


(A 16) 
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since 



Winged-Elliptic Cross Section with Semiminor Axis and Wing Planform 
Perpendicular to Crossflow Velocity 


For a winged-elliptic cross section 
where the semiminor axis b and the wing 
planform are perpendicular to the cross- 
flow velocity V n (sketch (o)), the same 
procedure is followed as for the previous 
configurations. 

For an elliptic body as shown in 
sketch (o), 





= 1 


(A 17) 


tan 5 = 


dy = 

dx 




and 


tan 2 6 _ b 2 - y 2 

1 + tan 2 5 b 2 - y 2 \ 1 - {a 2 lb 2 )] 


(A18) 


(A19) 


With equation (A19) combined with equations (Al), (A2), and (A7), we obtain 




4 f b b 2 - y 2 
dj b 2 - y 2 [\ - (a 2 /b 2 )] 



(A20) 


where d = l^fab is the diameter of the equivalent circular cross section. As before, the cross- 

flow drag coefficient for the equivalent circular cross section is given by Cj = 4/3. 

a n 

Thus, by integrating equation (A20), substituting 2 y/ah for d, and dividing the result 
by 4/3, we obtain 
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APPENDIX B 


DERIVATION OF (CJC„ ) FOR WINGED-SQUARE CROSS SECTIONS 

° Newt 

WITH ROUNDED CORNERS 


For winged-square cross sections with 
rounded corners (sketch (p)). the pressure coef- 
ficient over the front face is given by 

Cp = 2 sin 2 5 (Bl) 


r = kw 


V n 



where 5 is the local angle that a tangent to a 
forward-facing surface makes with the free- 
stream direction. In Newtonian flow, the pres- 
sure coefficients over the rearward face are 
assumed to be zero. The total section drag coef- 
ficient (based on width w) is then given by 



Sketch (p) 



For a rounded corner, as shown in sketch (p), 


(B2) 


and 


.v 2 -f y 2 = r 2 


tan 5 = 


dy 

dv 



sin 2 5 = 


tan 2 5 
1 + tan 2 5 



(B3) 

(B4) 


(B5) 
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The pressure coefficient for the rounded corner, obtained by substituting equation (B5) into 
(Bl), is 


C r 


corner 



(B6) 


The pressure coefficients for the front flat portion of the body and for the wing are 
given by 


Cp — 2 sin 2 5 = 2 

since 6 = 90°. 

Thus, with equations (B6) and (B7) substituted into equation (B2), we write 


C A = 


w 


/ 


(w/2 )~r 


dy + 


YV 


o o 

From integration of equation (B8), we obtain 

4 


f / jr 

J o w/2 




(B7) 


(B8) 


(B9) 


Now, to obtain (C„/C„ ) , it is necessary to find the diameter d of the equivalent 

n o Newt 

circular cross section which has the same area as the cross section studied. The area of the 
equivalent circular cross section is, of course, 

A = 4 d 2 (BIO) 

4 

and the area of the cross section studied is 

A = w 2 [ 1 - k 2 (4 - 7 r ) ] (B I 1 ) 


By equating equations (B 1 0) and (B 1 1 ), we obtain the equivalent diameter 

d = 2w 




it)k 2 


With C j from equation (B9) based on d instead of w, 

(M) 




k\ w 
d 


(Bl 2) 


(B 1 3 ) 
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where, for the equivalent circular section, 



(B14) 


Thus, from equations (B 1 2), (B13), and (B14), we obtain 




3) ]/'l - (4 - 7T)* 2 : 


0 < k < 0.5 


(B 1 5) 
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APPENDIX C 


FORMULAS TO COMPUTE GEOMETRIC PARAMETERS FOR TANGENT OGIVES 


To compute the aerodynamic characteristics of bodies of revolution having tangent 
ogive nose shapes, various geometric parameters must be obtained. To compute Cyy and C m , 
the planform area A p is required; to compute C m , it is also necessary to obtain volume V 
and distance x from the nose vertex to the centroid of planform area. To compute skin- 
friction drag, the wetted surface area A s is needed. 

For an ogival nose of length Cjy and diameter d (see sketch (q)) the following useful 
formulas have been derived (ref. 9): 

) 

(Cl) 


(C2) 


Sketch (q) 



_ i N (2/3) |/? 3 - [R* - {SL N ld?f /2 \ - (& N /d) 2 [ R - (1/2)] 
d ~ d A p /d 2 



where R is the ratio of the ogival arc radius r a to base diameter d and 



(C3) 


(C4) 


(C5) 
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Variation of crossflow drag coefficient with crossflow Mach number for circular cylinders 






Figure 2.— Variation of crossflow drag coefficient with crossflow Reynolds number for circular cylinders 

at subcriticai crossflow Mach numbers {M n < 0.4). 







Figure 3.- Variation of crossflow drag coefficient with crossflow Reynolds number for 
circular cylinders at supercritical Reynolds numbers and at crossflow Mach numbers from 
0.25 to 0.50 (from ref. 54). 



Circular cylinder length-to- diameter ratio 
Flat plate length - to -width ratio 

Figure 4.- Ratio of crossflow drag coefficient for a finite-length cylinder (or flat plate) to 
that for an infinite-length cylinder (or flat plate) (from ref. 56). 
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Figure 5. Variation of p ( ( / with M f[ (crossflow Mach number) computed from expert 
mental C y data (ref. 1 6) tor bodies of revolution at high a (45° < a < 60°). 
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Figuie 6. Variation of p with obtained Irom experimental results 

in figures 1, 4, and 5. 
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(a) Semimajor axis a perpendicular to crossflow velocity V n . 
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(b) Semiminor axis b perpendicular to crossflow velocity V fJ . 

. Ratio of local normal-force coefficient for an elliptic cross section 
to that for the equivalent circular cross section. 










Figure 8.— Ratio ofloeal normal-force coefficient for a winged-circular cross section 
to that for the equivalent circular cross section. 
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Figure 9.- Ratio of local normal-force coefficient fora winged-elliptic cross section to that 
for the equivalent circular cross section; semimajor axis a perpendicular to crossflow 
velocity V n . 
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Figure 10.- Ratio of local normal-force coefficient fora v 
for the equivalent circular cross section: semiminor ; 
velocity V . 
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Figure 1 2.- Bodies for which the aerodynamic characteristics were measured in reference 63 and 

computed in the present study for M 00 = 2.9. 



Computed Measured 

( ref. 63} 



Figure 13.— Comparison of computed with measured aerodynamic characteristics for cylindrical 
bodies with conical noses of fineness ratio 3 ; = 2.9. 




bodies with ocival noses of fineness ratio 5; M { 





Computed Measured 



Figure 16.- Comparison of computed with measured aerodynamic characteristics for cylindrical 
bodies with conical and ogival noses of fineness ratio 3; = 2.9. 








(a) Cylindrical bodies with conical 
noses of fineness ratio 3. 




(b) Cylindrical bodies with ogival noses 
of fineness ratio 5. 


(c) Cylindrical bodies with conical and 
ogival noses, overall fineness ratio 9. 


Figure 1 7.- Comparison of computed with measured aerodynamic force centers for cylindrical 
bodies with conical and ogival noses; M = 2.9. 
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Figure 18.- Ogive-cylinder bodies for which the aerodynamic characteristics have 
been computed and measured for 0.6 < M < 2.0. 


89 






iZ 


90 


gure 19 - Effect of Mach number on comparison of computed with measured aerodynamic 
characteristics for body V-C, : 0.6 0 < 2.0. 
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gure 20.“ Effect of Mach number on comparison of computed with measured aerodynamic 
characteristics for body A 7 C \ ; 0.6 <A/ 0 0 < 2.0. 





92 


Figure 21 Effect of Mach number on comparison of computed with measured aerodynamic characteristics 

for body N 2 C, C 2 ; 0.6 < A/,*, < 2.0, Re = 4.3X 1 0 5 . 






Figure 21.- Continued 




Figure 21.- Continued 







Figure 21 Concluded 
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Figure 22.- Effect of Mach number on comparison of computed with measured aerodynamic characteristics 

forbodyiV 7 C,C 3 ;0.6<iW oo <2.0, Re = 4.3X10 5 . 




Figure 22.- Continued 
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Figure 22. Concluded. 
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PLANFORM AND END VIEWS OF BODIES 
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Figure 24.— Bodies of elliptic cross section for which the aerodynamic characteristics were measured 

in reference 18 and computed in this study. 
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'O n f igurations studied. 
Concluded. 
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(b) Planform view at 0 - 90°. 



(c) Three-quarter view. 

Figure 25.- Photographs of aluminum body B 3 . 
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CONFIGURATION COMPUTED MEASURED 
PLAN VIEW 
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CONFIGURATION COMPUTED MEASURED 
PLAN VIEW 
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Figure 26.— Continued. 





CONFIGURATION COMPUTED MEASURED 
PLAN VIEW 





igure 26.— Continued. 
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PLAN VIEW 
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CONFIGURATION COMPUTED MEASURED 
PLAN VIEW 
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Figure 27.— Comparison of computed with measured aerodynamic characteristics for bodies 

and B 3 ; 0.6 < <2.0. 
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PLAN VIEW 



-igure 27.— Continued. 
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CONFIGURATION COMPUTED MEASURED 
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PLAN VIEW. 
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Figure 27.- Concluded. 
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Figure 28.- Comparison of computed with measured aerodynamic characteristics for body 

B 3 at 0 = 90°: 0.6 <^<2.0. 





Figure 28.— Continued. 
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Figure 28.- Continued. 





Figure 28.- Continued. 





Figure 28. - Concluded. 














(c) Body B 2 with wings W 2 , W A , and W s . 
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Figure 29.-- Concluded. 





Figure 30.— Plantorm views ot configurations for which the aerodynamic characteristics 
were measured in references 20 and 21 and computed in this study. 
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Figure 31.— Comparison of computed with measured aerodynamic characteristics for 

body-wing model R x W 5 . 
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Figure 3 1 Continued. 
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Figure 31.— Continued. 
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Figure 31.— Concluded 
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Figure 32.- Comparison of computed with measured aerodynamic characteristics for 

body-wing model B i W 2 . 
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Figure 32.— Concluded. 
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Figure 33.— Comparison oi computed with measured aerodynamic characteristics for 

body-wing model B 2 W s . 
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Figure 33.— Concluded. 
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Figure 34. Comparison of computed with measured aerodynamic characteristics for 

body-wing-tail model B } W 5 T. 
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Figure 34.— Concluded. 


CONFIGURATION 

PLAN VIEW METHOD COMPUTED MEASURED 
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Figure 35.— Comparison of computed with measured aerodynamic characteristics for 

body-wing-tail model R x W 2 T. 
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Figure 35.“- Concluded. 
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Figure 36.- Comparison of computed with measured aerodynamic characteristics for 

body-wing-tail model B 2 W 5 T. 
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Figure 36.— Concluded. 






0.5 d h -- l M = 3 d - 



Figure 37.- Additional tangent ogive noses and strake for modification of body B } . 



Figure 38.- Plan-view sketches of configurations used for vapor-screen and 

oil-flow studies. 
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6 FOOT W.T. VAPOR SCREEN APPARATUS 
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(a) Schematic drawing of equipment. 



(b) Picture of light-box installation. 


Figure 39.— Continued. 
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(c) Picture of model and camera installation. 


Figure 39.- Concluded. 
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STATION I 
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Figure 40.- Identification of typical items shown in vapor-screen photographs taken for a 

body at = 2.0. 
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Figure 41.— Vapor-screen photographs for B x C\ (body with fineness-ratio-3 nose). 
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Figure 41.— Continued. 





STATION 
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Figure 42.— Vapor-screen photographs for N 2 C\ (body with fineness-ratio-3.5 nose). 
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Figure 42.- Continued. 
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Figure 43. Vapor-screen photographs for N 3 C X (body with fineness-ratio-5 nose) 



tinued. 
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Figure 44.— Vapor-screen photographs for (V 4 Ci (body with blunt nose of fineness ratio 3). 
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Figure 45.- Vapor-screen photographs for A'sCjS (body with fineness-ratio-5 nose and afterbody strake). 
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Figure 45.— Continued. 



oncluded. 
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Figure 46.- Continued. 


STATION I 



Figure 46.— Concluded 
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Figure 47.- Vapor-screen photographs for B 2 at (p = 90° (body with constant a/b cross sections). 
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Figure 47.— Concluded. 
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Figure 48.— Continue. 




t 0 = 90° (body with variable a/b cross sections). 


STATION I 



Figure 49.— Continue 
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Figure 50. — Vapor-screen photographs for B X W S = N x C x W 5 (body with fineness-ratio-3 nose and aspect-ratio-3 wing) 
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Figure 5 1 Continued. 
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Figure 5 1 Concluded. 
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Figure 52.— Vapor-screen photographs for B { W 2 T = N x C x W 2 T (body with fineness-ratio-3 nose, aspect-ratio^ wing, and tail). 
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Figure 52.— Concluded 
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Figure 53.- Continued. 


STATION I 



Figure 53,— Concluded 
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Figure 54.- Vapor-screen photographs for W 2 T (body with fineness-ratio-5 nose, aspect-ratio-4 wing, and tail). 
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Figure 54.- Continued. 


STATION I 



Figure 54.— Concluded 
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(a) = 0.6. 

Figure 55.- Oil-flow photographs for/tj = A',C , (body with fineness-ratio-3 nose). 
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(b) = 0.9. 

Figure 55. Continued. 


TOP 


SIDE 


183 




(c) A/oo =2.0. 
Figure 55.- Concluded. 
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a «40° 



(a) M oo = 0.6. 

Figure 56. Oil-tlow photographs for N 2 C } (body with fineness-ratio-3.5 nose). 
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(b) = 0.9. 

Figure 56.- Continued. 
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TOP 






(c) M 0 o = 2.0. 
Figure 56.- Concluded. 
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(a) A/qq - 0.6. 

Figure 57.- Oil-flow photographs for /V 3 C, (body with fineness-ratio-5 nose). 
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TOP 


SIDE 





(b) A*,*, = 0.9. 
Figure 57.- Continued. 
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a *40° 




(c) Moo ~ 2.0. 
Figure 57.- Concluded. 
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a «40° 



(a) Mqo = 0.6. 

Figure 58.- Oil-flow photographs for N 4 C X (body with blunt nose of fineness ratio 3). 
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at » 30° 





(b) = 0.9. 

Figure 58.— Continued. 
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TOP 






(c) A/qo = 2.0. 
Figure 58. Concluded. 
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(a) Mqo = 0.6. 

Figure 59. Oil-How photographs for JV 3 C’,S (body with fineness-ratio-5 nose and afterbody strake). 
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(b) *00 = 0.9. 
Figure 59.— Continued. 
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(c) M 00 = 2.0. 
Figure 59.— Concluded. 
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(a) A/qq = 0.6. 

Figure 60.— Oil-How photographs for B 2 at <p = 0° (body with constant alb cross sections). 
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(c) A/qo = 2.0. 
Figure 60.- Concluded. 
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(a) /Woo = 0.6. 

Figure 61 .- Oil-flow photographs for B 2 at 0 = 90° (body with constant a/b cross sections). 
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(b) A/,*, =0.9. 
Figure 61.- Continued. 
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a«30° 




(e) Moq - 2.0. 

Figure 61. Concluded. 
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(b) Al 00 = 0.9. 
Figure 62. Continued. 
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SIDE 


a«IO° 






(c) M 00 = 2.0. 
Figure 62. Concluded. 
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(a) M ^ = 0.6. 

Figure 63. Oil-How photographs for B 3 at </> = 90° (body with variable a/b cross sections). 
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(b) /Woo — 0.9. 
Figure 63. Continued. 
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(c) Moo = 2 . 0 . 
Figure 63. — Concluded. 
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(a) M „ = 0.6. 

Figure 64. Oil-How photographs for W 2 T — (V, C { W 2 T (body with fineness-ratio-3 nose, 

aspect-ratio-4 wing, and tail). 
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POST-RUN SAME AS 
PRE-RUN OIL PATTERN 


(a) A/qo = 0.6 — Concluded. 
Figure 64. - Continued. 
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(b) Mao ~ 0-9. 
Figure 64. - Continued. 
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-POST-RUN SAME AS 
PRE-RUN OIL PATTERN 
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POST-RUN SAME AS 
PRE-RUN OIL PATTERN 


(b) Mqq = 0.9 - Concluded. 
Figure 64.- Continued. 
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(c) M oa = 2.0. 
Figure 64. Continued. 
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PRE-RUN OIL PATTERN 


(e) A/qq = 2.0 — Concluded. 
Figure 64. Concluded. 
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(a) Moo = 0.6. 


Figure 65. Oil-tlow photographs for A’ 3 C, W 2 T (body with fineness-ratio-5 nose, 

aspect-ratio-4 wing, and tail). 
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(a) M x = 0.6 Concluded. 


Figure 65.— Continued. 
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TOP 


SIDE 



(b) M oe> = 0.9. 
Figure 65.- Continued. 
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a » 30° 
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a « 40° 



-7 POST- RUN 
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OIL PATTERN 


(b) M 00 = 0.9 - Concluded. 


Figure 65. Continued. 
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(c) M ao = 2.0. 
Figure 65.-- Continued. 
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POST -RUN 
SAME AS 
PRE-RUN 
OIL PATTERN 


(c) ^00 = 2.0 Concluded. 


Figure 65.— Concluded. 
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ogive-cylinder bodies of revolution. 
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Figure 66.— Continued 
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Figure 67.— Concluded. 
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for an ogive-cylinder body of revolution. 
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Figure 69. Effect of afterbody side strakes on side-force coefficient and side-force position for an 
ogive-cylinder body of revolution with a fineness-ratio-5 nose. 
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Figure 69.— Concluded. 
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Figure 70.- Effect of elliptic cross section and roll angle on side-force coefficient and 
side-force position (elliptic body of constant a/b). 
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Figure 70.- Concluded. 
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Figure 7 1 Effect of elliptic cross section and roll angle on side-force coefficient and 
side-force position (elliptic body of variable a/b). 
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Figure 7 1 Concluded. 

















